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ABSTRACT 

Performance  data  from  the  NASA-Army  OH-6A  higher  harmonic  control 
(HHC)  flight  test  program  showed  significant  reductions  in  main  rotor  shaft  torque 
and  engine  power  in  hover  and  forward  flight  [1].  The  unsteady  aerodynamics  with 
the  higher  harmonic  control  system  application,  including  wake  effects,  were  con- 
sidered   to  study  whether  such  power  reductions  are  feasible. 

An  airfoil  oscillating  in  pure  plunge  can  achieve  propulsive  force  ("Katzmayr 
effect"),  as  in  the  case  of  birds'  wings  flapping  in  tlight.  Here  it  will  be  shown 
that   this  effect  can  be  enhanced    in  the  presence  of  layers   of  shed   vorticity  with 
the  proper  phasing.     In    addition,  while  it  is  known    that  an  airfoil  oscillating  in 
pitch,  will  typically  produce  drag  at  most  values  of  reduced  frequency,  it  is  found 
that  the  presence  of  another  layer  of  shed    vorticity    of  the  proper  phase,  can 
reduce  the  drag  on  the  pitching  airfoil  depending  upon  wake  spacing,  reduced  fre- 
quency, and  phase.    Under  some  conditions  the  added  layer  or  layers  of  shed  vor- 
ticity   will  even  result  in  propulsive  force  acting  on  the  pitching  airfoil  similar  to 
the  "Katzmayr"  effect  for    the  plunging  case. 

It  was  found,  for  the  OH-6A  helicopter,  that  the  measured    reductions  in 
main  rotor  shaft  torque  and  engine  power  are  feasible  when  evaluated  with   respect 
to  the  "Katzmayr  effect"  and  the  additional  drag  reduction  or  propulsive  force 
obtained    due  to  pitch    and    plunge  oscillations  with  the  effect  of   adjacent  wake 
layers  of  shed  vorticity. 
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Summary 

Open  loop  performance  data  from  the  NASA-Army  OH-6A  higher  harmonic 
control  (HHC)  flight  test  program.  Wood  et  al.  [1],  show  significant  reductions 
in  main  rotor  shaft  torque  and  engine  power  in  the  airspeed  regime  from  hover 
to  100  knots.  Depending  upon  HHC  controller  phase  and  helicopter  airspeed, 
reductions  in  power  were  recorded  as  large  as  20%,  with  reductions  of  the  order 
of  10  %  being  typical. 

Higher  harmonic  control  is  an  active  helicopter  vibration  control  concept 
which  alters  the  aerodynamic  loads  on  the  rotor  blades  such  that  the  blade  re- 
sponse is  reduced.  This  in  turn  reduces  the  vibratory  forces  and  moments  acting 
at  the  hub,  which  cause  airframe  vibration.    Basically,  HHC  is  an  electronic, 
computer-controlled  active  vibration  suppression  system  which  senses  and  cancels 
vibrations  in  a  helicopter  by  N/revolution  feathering  or  pitch  motion  of  the  rotor 
blades,  N  being  the  number  of  blades  .     Very  substantial  vibration  reduction  was 
achieved  under  a   NASA-Army  sponsored  program  using  a  modified  OH-6A  he- 
licopter. Higher  harmonic  blade  pitch  control  was  implemented  by  superimposing 
4/rev.  (32  Hz)  swashplate  motion  upon  basic  cyclic  and  collective  control  inputs. 


K.Roberts  Wood.  R.  W    Powers.  J  H.  Cline.  and  C  E.  Hammond.    "On    Developing  and 
Flight  Testing  a  Higher  Harmonic  Control  System."   Journal    ofA.H.S  January.    1085. 
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During  the  program,  in  addition  to  reducing  vibration  ieveis,  the  HHC  system 
showed  the  potential  for  decreased  helicopter  power  requirements.  Recent  study 
of  the  OH-6A  flight  test  data  indicates  that  the  mechanism  by  which  the  power 
reduction  was  achieved  is  related  to  the  unsteady  aerodynamics  associated  with 
HHC.  This  is  supported  by  the  fact  that  higher  harmonic  control  by  the  very 
nature  of  the  method  by  which  it  is  implemented  through  oscillating  an  airfoil 
in  pitch  with  resulting  plunge  motion,  requires  definition  of  unsteady  aerodynam- 
ics to  adequately  model  its  physics 

The  fundamental  closed  form  solutions  of  Theodorsen  [3],  and  Loewy  [4] 
provide  the  basis  for  theoretical  work    in  this    area.     The  closed    form  theory 
shows  rapid  changes  in  the  lift  deficiency  function  with  changes  in  reduced  fre- 
quency, wake  spacing  and  frequency  ratio.    In  the  past,  emphasis  in  the  study 
of  unsteady  aerodynamics  has  tended  to  focus  on  flutter  instability  and  the  effect 
of  unsteady  aerodynamics  on  generating  lift  and  torsional  loads. 

In  this  research,  we  are  interested  in  performance,  and  emphasis  will  be  on 
the  effect  of  unsteady  aerodynamics  on  the  drag  of  the  airfoil.  When  the  drag 
force  is  reversed,  acting  to  propel  the  airfoil  forward,  it  is  generally  referred  to 
as  either  negative  drag  or  propulsive  force.  The  classic  reference  on  this  subject 
is  that  by  Garrick  [5] 


1   E.R.Wood.    Max  F.  Platzer.    Ahmed  Abourahma.    Mark    A.  Couch  "  On  the  Unsteady 
Aerodynamics    of   Higher  Harmonic  Control."  Paper  No    C  17.  Nineteenth  European 
Rotorcraft  Forum.  Cemobbio  (Como).  Italv.   September  1993. 
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While  Garrick's  [5]  work  shows  that  an  airfoil  oscillating  in  pitch  will 
typically  produce  drag  in  the  lower   reduced  frequency  range,  k.  it  is  found  that 
the  presence  of  another  layer  of  shed  vorticity  of  the  proper  phase  can  reduce 
the  drag  on  pitching  airfoils  and  depending  upon  wake  spacing,  reduced  frequen- 
cy, and  phase,  may  even  enhance  the  propulsive  force  acting  on  the  airfoil  similar 
to  the  "Katzmayr  "  effect    for  plunging  airfoils. 

The  specific  mechanisms,  that  have  been  used    in  this  research  to  discuss 
me  phenomenon  of  OH-6A  helicopter  rotor  power  reduction,  are:  (  1 )  the  unsteady 
.^dynamics  associated  with  a  plunging  airfoil  and  how  that  produces  a  propul- 
sive force  or  decrease  in  rotor  power;  (2)    the  influence  of  a  layer  of  wake 
vorticity  on  this  propulsive  force,  which  represents  the  influence  of  the  unsteady 
aerodynamics  of  one  rotor  blade  on  adjacent  rotor  blades,  (3)  the  unsteady  aero- 
dynamics associated  with  an  airfoil  oscillating  in  pure  pitch  and  how  this  pro- 
duces in  general  a  drag  force  on  the  airfoil  unless  accompanied  by  layers  of  shed 
wake  vorticity ;  and  (4)  the  influence  of  the  phasing  of  a  layer  of  shed  vorticity 
on  this  drag  force,  and  how,  depending  upon  specific  phasing,  the  effect  on  the 
airfoil  will  be  either  drag  or  propulsive  force. 

The  analysis  of  the  above  effects  was  based  upon;  (1)  the  classical  flat 
plate  unsteady  aerodynamic  theory  for  an  incompressible  fluid  of  Theodorsen  [3], 
its  extension  to  include  chordwise  forces  (drag)  by  Garrick  [5],    the  rotor  wake 
analysis  of  Loewy  [4],  (2)  the  panel  method  developed  by    Platzer  et  al.  [6], 
that  can  be  applied  for  unsteady  incompressible  flow  past  airfoils  or  airfoil  com- 
binations of  arbitrary  geometry  and  amplitudes  of  different  kinds  of  motions. 
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In  the  first  approach.  Theodorsen's  well  known  lift  deficiency  function  C(k) 
has  been  modified  by  Loewy  for  the  case  of  a  hovering  rotor.  Airfoil  drag  or 
propulsive  force  was  determined  by  the  methods  of  Garrick  [2],  In  the  second 
approach,  the  panel  code  allowed  a  systematic  study  of  the  lift  and  drag  produced 
by  oscillating  airfoils  and  airfoil  combinations  in  an  incompressible  inviscid 
flow. 

The  panel  code  results  were  checked  against  different  references.  The  results 
showed  good  agreement  with  these  references.  For  the  parameters  ot  the  OH- 
6A  rotor  in  hover,  the  main  conclusion  was  that  the  reductions  in  drag  and 
propulsive  forces  generated  on  an  airfoil,  oscillating  in  pitch  or  plunge,  due  to 
wake  phasing,  are  sufficient  to  explain  the  reductions  in  power  recorded  on  the 
OH-6A  during  the  NASA- Army  flight  test  program  [1], 
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I .  INTRODUCTION 

A.  GENERAL 

The  determination  of  the  unsteady  aerodynamics  associated  with  a  helicopter  rotor  has 
presented  a  challenge  to  the  helicopter  analyst  since  the  conception  of  the  helicopter.  The 
effect  of  blade  oscillations  on  resulting  airloads  has  prompted  new  interest  in  oscillatory 
blade  aerodynamics  in  connection  with  possible  performance  gains.  Rotarv  wing  unsteady 
aerodynamics  is  considerably  more  complex  than    fixed  wing  aerodynamics. 

The  simplest  case  is  that  of  hover,  but  even  here  the  problem  is  complicated  by  the 
influence  of  the  shed  helicopter  wake  from  the  rotor.  In  forward  flight,  the  problem  is 
further  complicated  by  the  fact  that  the  velocity  of  a  local  blade  element  ranges  at  the 
blade  tip  from  transonic  on  the  advancing  side  to  low  subsonic  on  the  retreating  side,  and 
at  the  blade  root  from  low  subsonic  on  the  advancing  side  to  reversed  flow  on  the  retreat- 
ing side. 

For  this  reason  there  is  a  great  need  for  advances  in  the  present  theory  in  order  to 
explain  modern  helicopter  problems,  yet,  the  number  of  advances  in  rotorcraft  unsteady 
aerodynamic  theory  are  very  limited,  so  we  need  to  review  current  unsteady  aerodynamic 
theory  for  rotorcraft  and  to  advance  the  theory  in  one  of  the  areas  where  improvement  is 
badly  needed. 

The  goal  of  this  thesis  is  to  explain  the  measured  power  reductions  in  the  OH-6A 
higher  harmonic  control  (HHC)  flight  test  program  conducted  by  E.  R.  Wood  et  al  [1]. 
This  research  will  consider  the  unsteady  aerodynamics  with  HHC  system  application, 
including  wake  effects,  in  order  to  determine  whether  such  power  reductions  are  feasible. 
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B.  ROTARY  WING  UNSTEADY  AERODYNAMICS 

In  this  section,  some  elements  of  unsteady  airfoil  aerodynamics  will  be  discussed. 
Consider  an  airfoii  undergoing  sinusoidal  pitch  and  plunge  motion.  As  the  airfoil  performs 
an  oscillation,  vortices  are  shed  into  the  medium  with  a  circulation  strength  equal  in  mag- 
nitude to  the  increase  in  circulation  about  the  airfoil,  but  opposite  in  direction.  These  dis- 
turbances are  stored  in  the  fluid  because  the  shed  vorticitv  convects  downstream  at  the 
local  flow  field  velocity. 

The  counter-rotating  vortices  induce  a  sinusoidal  tlow  field  which  further  changes  the 
X  The  result  of  this  complex  flow  field  is  a  time  difference  or  delay  between  the 
airfoil's  motion  and  the  induced  aerodynamic  forces.  This  delay  is  known  as  the  phase 
lag,  y\j  [34]  .  To  simplify  calculations  for  this  type  of  motion,  it  is  common  to  describe 
the  airfoil  position  and  the  associated  aerodynamic  forces  by  complex  variables.  For  pure 
plunge  oscillations  the  vertical  motion  of  the  airfoil  is  described  by  the  real  part  of  the 
following  equation: 

h(t)   =   h0    elcot    (1.1) 

where  h0    is  a  complex  number,  and  o  is  the  frequency  of  oscillation.  Similarly,  the  lift  is 
dec~      d  by  : 


L=L0reitot   (1.2) 


where  L0    is  the  quasi-steady    lift  given  by  the  expression 


L0=.5  pU2S  C!a  (h  /  U)  (1.3) 
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This  is  termed  quasi-steady  because  the  angle  of  attack  is  represented  by  h/U. 
The  values  of  r0  and  \\i  represent  the  magnitude  and  the  phase  respectively,  of  the 
true  instantaneous  lift  relative  to  the  quasi-steady  lift.  The  variables  r0  and  \u  in 
general  depend  on  the  reduced  frequency  k,  the  Mach  number  M.  and  the  Reynolds 
number.  For  inviscid.  incompressible  fluid  the  values  of  ru  and  \u  will  be  only  a 
function  of  k  [34],  A  complete  solution  for  the  oscillating  flat-plate  airfoil  in  an 
incompressible  inviscid  flow  has  been  obtained  by  Kussner  and  Theodorsen  and 
reproduced  from  Reference  [35]  in  Figure  (1.1  ). 


Figure  (1.1)      l^/l-,,  Wiase  Diagram  (   Imaginary     vs   Real  ) 


C.  ADDITION  OF  WAKE   EFFECT 

For  the  helicopter  in  hover,  the  rotation  of  the  wing  introduces  a  number  of  features 
that  require  special  attention,  notably  the  returning  vortex  wake,  the  time  varying  free 
stream,  and  radial  flow,  and  a  fundamentally  transcendental  geometry  that  requires  either 
approximate  or  numerical  solutions.  The  lift  on  a  wing  is  due  to  its  bound  circulation. 
Conservation  of  vorticity  in  the  tlow  requires  that  there  be  trailed  and  shed  vorticity  in 
the  wake  behind  the  wing. 

The  spanwise  variation  of  bound  circulation  results  in  trailed  vorticity  parallel  to  the 
free  stream  direction.  Time  variation  of  the  bound  circulation  leads  to  shed  vorticity  par- 
allel to  the  wing  span  The  wake  is  composed  of  sheets  of  vorticity  convected  downstream 
from  the  trailing  edge  by  the  free  stream  velocity.  The  wake  of  a  rotor  in  hover  or  in  ver- 
tical flight  consists  of  helical  vortex  sheets  below  the  disk,  one  from  each  blade.  Unsteady 
motion  of  the  rotor  blade  will  produce  a  shed  vorticity  in  the  wake  spirals.  With  low  disk 
loading  ( at  low  collective  pitch  settings  ),  the  wake  remains  near  the  rotor  disk  and  there- 
fore passes  close  to  the  following  blades.  Thus  the  wake  vorticity  is  not  convected  down- 
stream of  the  airfoil  as  with  a  fixed  wing,  and  the  shed  vorticity  sheets  below  the  rotor 
disk  must  be  accounted  for  to  correctly  estimate  the  unsteady  loads. 

For  high  inflow  (high  collective  pitch  settings  )  or  in  high-speed  forward  flight,  the 
wake  will  be  convected  away  from  the  blades.  The  returning  shed  wake  influence  is  pri- 
marily a  concern  of  hover  and  vertical  flight,  and  forward  flight  at  airspeeds  less  than 
100  knots.  Assuming  high  aspect  ratio  of  the  blade,  lifting  line  theory  requires  a  knowl- 
edge of  the  loads  on  the  blade  section,  and  the  returning  shed  wake  of  the  rotor  must  be 
incorporated  into  the  two-dimensional  unsteady  airfoil  theory. 


The  wake  far  from  the  blade  section  will  have  little  influence,  so  the  emphasis  is  on 
modeling  the  wake  near  the  blade,  which  for  low  inflow  consists  of  vortex  sheets  that  are 
nearly  planar  surfaces  parallel  to  the  disk  plane.  Based  on  these  considerations,  a  two- 
dimensional  model  for  the  unsteady  aerodynamics  of  the  rotor  can  be  constructed.  Loewy 
in  1957,  developed  a  two-dimensionai  model  for  the  unsteady  aerodynamics  of  the  blade 
in  a  hovering  rotor.,  including  the  effect  of  the  helical    shed  wake. 

D.    HIGHER  HARMONIC  CONTROL  (HHC) 

The  primary  source  of  helicopter  vibration  is  the  higher  harmonic  blade  loads  gener- 
ated by  the  aerodynamic  environment  at  the  rotor  disk.  For  an  N-bladed  rotor,  the  oscil- 
latory (N-l)P,  NP,  and  (N+1)P  blade  loads  are  transmitted  as  exciting  hub  forces  and 
moments  to  the  airframe  at  a  dominant  frequency  of  N  per  rev  (NP).  A  significant  reduc- 
tion in  vibration  level  has  been  achieved  by  industry  in  the  last  twenty  five  years.  The 
desired  goal  was  to  achieve  0.02  g  and  this  goal  cannot  be  met  without  a  quantum 
advance  in  vibration  control  technology.  The  vibration  levels  first  specified  by  the  U.  S 
.Army  for  the  AAH/UTTAS  procurement  in  the  early  1970s  were  later  revised  upward  in 
the  mid  1970s  to  better  reflect  realistic  design  goals  consistent  with  the  state-of-the-art  in 
helicopter  vibration  control  [1]. 

A  higher  harmonic  control  system  drives  the  blades  in  pitch  at  the  (N-l)P,  NP,  and 
(N+1)P  harmonics  of  rotor  rotational  speed,  generating  new  unsteady  airloads,  which  in 
combination  with  the  new  oscillatory  inertial  loads,  cancel  the  harmonics  of  blade  loads 
that  cause  airframe  vibration.  Therefore  the  vibrations  are  suppressed  at  the  source. 
Although  there  are  several  ways  to  implement  HHC  on  a  rotor  system,  the  approach 
generally  followed  to  date  has  been  by  blade  root  feathering  using  swashplate  oscillations . 


By  means  of  electro-hydraulic  servo-actuators,  the  svvashplate  is  excited  in  the 
collective,  longitudinal  cyclic,  and  lateral  cyclic  modes  at  NP  resulting  in  blade  pitch 
oscillations  at  three  distinct  frequencies  of  (N-l)P,  NP,  and  (N+1)P  in  the  rotating  frame. 
In  the  system  considered,  developed  for  the  four-bladed  OH-6A  helicopter,  higher  har- 
monic blade  feathering  for  vibration  reduction  is  achieved  by  superimposing  4/rev  swash- 
plate  motion  upon  basic  collective  and  cyclic  flight  control  inputs. 

Perturbing  the  stationary  svvashplate  at  4/rev  both  vertically  and  in  pitch  and  roll 
results  in  third,  fourth,  and  fifth  harmonic  blade  feathering  in  the  rotating  system.  Fourth 
harmonic  blade  feathering  is  achieved  by  oscillating  the  svvashplate  vertically  about  its 
collective  position,  while  third  and  fifth  harmonic  blade  feathering  in  the  rotating  system 
results  from  4/rev  tilting  of  the  stationary  swashplate  in  pitch  and  roll  about  its  cyclic  tilt 
position  [2].  Results  of  experimental  efforts  by  Hammond  [52]  showed  that  successful 
suppression  of  vibration  can  be  achieved  by  oscillating  the  blades  at  relatively  small 
angles. 

E.  OH-6A  (NASA/ARMY/HUGHES)  TEST  PROGRAM  DATA 

The  goal  of  this  research  has  been  a  systematic  study  of  the  OH-6A  helicopter 
HHC  flight  test  program,  in  order  to  understand  and  explain  the  significant  reduction  in 
main  rotor  shaft  torque  and  engine  power  that  was  measured  in  flight.  Shown  in  Figure 
(  1.2  )  are  the  power  reduction  results  from  the  HHC  Open  Loop  Flight  Test  Program 
(NASA/  Army/  Hughes-  Sept.    1982). 
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Figure  (  1.2  ):  Power  Reduction  Results  from   Higher  Harmonic  Control    (HHC) 
Open  Loop  lest  Program  (NASA/Army/I  I  ughcs  -  Sept.  1982) 
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Figures  (13)  and  (1.4)  (From  Ref.  [1])  show  the  very  substantial  vibration  reduction 
achieved  under  a  NASA-Army  sponsored  program  using  a  modified  OH-6A  helicopter.  Higher 
harmonic  blade  pitch  control  was  achieved  by  superimposing  4/rev.  (32  \\z)  swashplate 
motion  upon  basic  evelic  and  collective  control  inputs.  The  aircraft  was  flown  from  zero 
airspeed  to  100  knots  with  the  HHC  system  operated  both  open  loop  (manually)  and  closed 
loop  (computer  controlled) 
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Figure  (  1.3  ):  Comparison  of  4/rev  Vertical  Vibration  at  the  Pilot 
Sent   Versus  Airspeed  with    HI1C    "on"  and   "off". 
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Figure  (  1.4  ):  Comparison  of  4/rev.   Lateral   Vibration  at  (lie  Pilot  Seat 
Versus  Airspeed  with   HFIC  "on"  and  "off. 
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The  significance  of  power  reduction  means  that:  ( 1 )  the  helicopter  payload  can  be 
increased.  (2)  the  maneuverability  can  be  markedly  increased.  (3)  marginal  operating  condi- 
tions, such  as  at  sea  in  high  sea  states  can  be  better  met.  and  (4)  the  helicopter  operating 
costs  can  be  substantially  reduced.  Increased  interest  in  HHC  and  its  potential  benefits  for 
the  helicopter  now  require  that  we  try  to  use  the  proper  tool  to  understand  this  phenome- 
non. It  would  appear  that  we  presently  have  the  computational  toois,  namely  the  panel 
method  for  unsteady  inviscid,  incompressible  flow  past  airfoils  or  airfoil  combinations  of 
arbitrary  geometry,  developed  by  Platzer  et  al  [6],  and  the  single  wake  closed  form  anal- 
ysis developed  from  Loewy's  theory  by  Wood  and  Couch  [40], 
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F.  SCOPE 

Chapter  II  contains  a  quick  review  of  the  most  important  research  starting  with  the 
work  of  Knoller[10]  in  1909  and  Betz[ll]  in  1912.  They  explained  the  bird's  ability  to 
generate  forward  thrust  by  means  of  wing  flapping.  In  1922.  R.  Katzmavr,  experimentally 
verified  the  work  of  Knoller  and  Betz.  Birnbaum,  von  Karman.  Theodorsen.  and  Garnck  in 
the  years   1924  through   1936.  were  the  pioneers    of   flapping    wing  aerodvnamics. 

The  fundamental  closed  form  solutions  by  Theodorsen  [3],  Loewy  [4],  Garnck  [5],  pro- 
vide our  basis  for  theoretical  work  in  this  field.  While  Loewy's  work  on  wake-induced 
flutter  helps  to  explain  the  phenomenon,  it  also  points  up  difficulties  to  be  overcome.  That 
is  the  closed  form  theory  shows  rapid  changes  in  the  lift  deficiency  function  with  changes 
in  reduced  frequency,  k;  wake  spacing,  h;  and  frequency  ratio,  m. 

Garrick  further  developed  the  work  of  Theodorsen  by  deriving  the  horizontal  force  for- 
mulas for  airfoils  or  airfoil-aileron  combinations  oscillating  in  any  of  the  three  degrees  of 
freedom:  vertical  flapping,  torsional  oscillation  about  a  fixed  axis,  and  angular  oscillation 
of  the  aileron  about  a  hinge.  Al so  presented  in  Chapter  II,  Loewy's  two-dimensional  model 
for  the  unsteady  aerodynamics  of  the  blade  of  a  hovering  rotor,  is  presented  to  provide  the 
base  for  discussing  the  obtained  results.  Also  presented  in  Chapter  II.  is  a  quick  survey  of 
the  most  relevant  references  who  developed  and  used  panel  methods  to  study  unsteady 
aerodynamics  of  oscillating  airfoils  in  an  incompressible,  inviscid  flow 
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In  the  third  chapter,  the  propulsive  force  is  estimated  according  to  the  OH-6A  helicop- 
ter parameters  obtained  from  reference  [36]  In  this  analysis,  the  blade  is  divided  into  sec- 
tions. Assuming  hovering  conditions,  the  reduced  frequency  at  each  section  was  evaluated 
Using  the  unsteady  panel  code  the  time  history  of  the  airfoil  drag  (  or  propulsive  force  ) 
coefficient  was  recorded.  For  the  assumed  plunging  motion  with  plunging  amplitude  as 
obtained  from  reference  [1],  for  a  given  mode  of  vibration,  the  average  propulsive  force 
for  each  section  was  calculated. 

The  totai  average  propulsive  force  was  obtained  by  adding  together  the  contribution 
from  each  blade  section.  The  obtained  results  showed  that  the  resulting  propulsive  force 
from  a  single  airfoil  oscillating  in  plunge  is  not  enough  to  explain  the  significant  power 
reduction  on  the  OH-6A  helicopter.  This  led  us  to  study  the  aerodynamic  interference 
between  two  blades,  which  is  the  main  subject  of  the  fourth  chapter. 

In  Chapter  IV,  the  two-airfoil  unsteady  panel  code  developed  by  Pang  [  25  ],  is  used 
to  study  the  wake  interference  between  two  airfoils.  This  code  has  been  formulated  to 
solve  for  the  potential  flow  for  two  airfoils  executing  unsteady  motion  in  an  inviscid 
incompressible  flow  medium.  It  is  an  extension  of  the  code,  developed  by  Teng  [24],  for 
single  airfoils.  The  technique  uses  the  well  known  panel  method  for  steady  flow  and 
extends  it  to  unsteady  flow  by  accounting  for  the  continuous  shedding  of  vortices  into  the 
trailing  wake. 
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Numerous  case-runs  are  presented  to  illustrate  the  aerodynamic  interference  between  two 
oscillating  airfoils.  Results  obtained  for  plunging  and  pitching  airfoils  show  a  significant 
production  of  propulsive  force  at  certain  wake  phase  angles.  In  Chapter  V  the  question  of 
propulsive  efficiency  due  to  plunging  or  pitching  and  the  possibility  of  single-degree-of- 
freedom  pitching  flutter  are  discussed.  Chapter  VI  discusses  the  findings  of  this  research 
and    presents  the  main  conclusions. 
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II.  BACKGROUND 

The  effect  of  plunge  and  pitch  oscillations  on  the  resultant  drag  of  an  airfoil 
has  been  the  object  of  aeroelastic  research  back  almost  to  the  time  of  the  Wright 
brothers  (1903)  In  1909  and  1912.  respectively,  Knoller  [10],  and  Betz  [11]  were 
the  first  ones  to  explain  the  birds"  ability  to  generate  a  forward  thrust  by  means 
of  wing  flapping  .  Another  one  of  the  earliest  investigators  of  this  era  was  R. 
Katzmayr  [7],  who  in  1922  conducted  two  series  of  experiments  to  explore  the 
effects  of  oscillatory  fluid  or  airfoil  motion  on  resultant  lift  and  drag. 

Katzmayr's  experiments  which  were  performed  in  the  aerodynamic  laboratory 
of  the   Technical  University  of  Vienna,  may  be  divided  into  two  series.  In  the  first 
series,  the  angle  of  attack  of  the  wing  model  was  changed  by  oscillating  the  wing 
model  about  an  axis  parallel  to  the  span  and  at  right  angles  to  the  airflow.  In  the 
second  series  of  experiments,  the  airfoils  were  stationary  while  the  airflow  itself 
was  subjected  to  periodic  oscillations.  Here,  for  certain  air  velocities  and  frequen- 
cies, and  considering    thick  airfoil  sections.  Katzmayr  measured  negative  drag  or 
positive  propulsive  force.  This  discovery,  to  be  later  known    as  "Katzmayr  effect", 
helped  to  explain  the  propulsive  mechanism  of  birds  flying  through  the  air. 


Platzer.  M.  F..  Neace.  K.  S..  and  Pang,  C.  K..." Aerodynamic  Analysis  of  Flapping  Wing 
Propulsion".    AIAA  paper  No.    93-0484.  31st    Aerospace  Sciences  meeting,  Reno.  NV,  January 
11-14.   1993. 
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Later.  E.  G.  Richardson  [8]  used  the  "Katzmayr  effect"  to  explain  the  locomotion 
of  fish  through  water.  In  his  paper,  Katzmayr  summarized  his  results  to  show:(l) 
the  effect  of  flowing  air,  whose  direction  is  undergoing  constant  periodical  change, 
is  quite  favorable,  (2)  the  wing  section  which  exhibits  favorable  characteristics  in  a 
constant  airflow,  works  better  in  an  oscillating  flow. 

In   1935,  Theodore  Theodorsen  [3]  introduced    the  general    theory    of   aerody- 
namic instability    and  the  mechanism  o^   flutter.   In    the  first  part  of   his    paper. 
Theodorsen  developed    the  velocity    potentials    due  to  flow  around    oscillating  air- 
foil-ailerons.   In  the  later  part  of  the  paper,  he  developed  the  differential  equations 
of  motions  and  determined  the  flutter  speed. 

The  most  relevant  work  with  respect  to  the  present  research  is  that  by  I.  E. 
Garrick  [5],  who  applied  the  work  of  Theodorsen  [3]  to  obtain  a  closed  form  ana- 
lytical expression  for  the  propulsive  force  generated  by  an  oscillating  airfoil.  An 
airfoil  oscillating  in  pure  plunge  or  flapping  motion  was  shown  by  Garrick  to  have 
a  propulsive  force  throughout  the  entire  range  of  reduced,  k,  while  for  an  airfoil 
oscillating  in  pitch,  propulsive  force  can  only  be  achieved  for  a  specified  range  of 
reduced  frequency,  k,  where  k  must  be  greater  than  2.0  . 
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In   1942,  W    Schmidt  [14]    recognized  that  the  Knoiler-Betz  or  Katzmayr  effect 
applies  to  both  a  flapping  airfoil  in  uniform  flow  or  a  stationary  airfoil  in  oscillat- 
ing flow.    He  demonstrated  experimentally  that  a  stationary  airfoil  positioned  in  the 
wake  of  a  flapping  airfoil  increases  the  propulsive  efficiency  to  almost  100  percent 
throughout  the  whole  frequency  range  because  the  stationary  airfoil  converts  the 
vortical  energy  generated  by    the  flapping  airfoil  into  additional  thrust.  He  also 
found  that  the  mechanical  limitations  inherent  in  pure  flapping  motions  could  be 
overcome  by  an  arrangement  which  he  called  the  wave  propeller. 

H.  Bosch  [15]  investigated  the  interference  of  two    lifting  surfaces  in  two- 
dimensional  unsteady    incompressible  flow,    using    classical  flat  -plate  theory.    He 
showed  that  a  harmonically  flapping  airfoil  upstream  of  a  stationary  airfoil  increases 
the  propulsive  efficiency  to  almost  100  %,  thus  confirming    Schmidt's  experimental 
findings.    Bosch's  analysis    is  limited  to  flat-plate  airfoils  oscillating  about  a  zero 
incidence  mean  position. 

In   1957,    Loewy  [4]  developed  a  two-dimensional  model  for  the  unsteady 
aerodynamics  of  the  blade  of  a  hovering  rotor.  He  derived  an  expression  for    the 
two-dimensional  oscillatory  lift  deficiency  function  as  a  function  of  the  reduced 
frequency,  wake  spacing  and  the  ratio  of  the  airfoil  frequency  of  oscillation  and  the 
rotational  speed. 
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The  wake  in  hover  or  vertical  flight  consists  of  helical  vortex  sheets  below  the 
disk,  one  from  each  blade.  Unsteady  motion  of  the  rotor  blade  will  produce  shed 
vorticity  in  the  wake  spirals. 

In    Figures  (2.1  A)  through   Figure  (2. ID),    from  Ref.  [4],    Loewy    illustrated  the 
basic  elements  of  the  unsteady    flow  field  of  a  helicopter  rotor.     A  rotor  in    hover- 
ing or    vertical  flight  produces  a  trailing    tip    vortex  with  a    downward  axial  veloc- 
ity   that,  if  otherwise  undisturbed,  forms  a  contracting  helix  as  shown  in  Figure 
(2.1  A).     If  the  inflow  over  the  rotor  disk  is   constant,  then  the  fluid   off  the  trailing 
edge  of  the  blades  makes  a  helical  surface  with  horizontal  radial  elements  as  shown 
in  Figure  (2.  IB). 

If  there  is  an  oscillation  in  blade  effective  angle  of  attack,  blade  lift  will  alter- 
nate also,  and  as  a  result  of  these  changes  in  lift,  vortices  will  be  shed  continu- 
ously at  the  blade  trailing  edge.  These  vortices  fall  along  the  horizontal  radial 
elements  of  the  helical  surface  shown  in  Figure  (2. IB),  so  long  as  the  oscillations 
in  angle  of  attack  are  small.    Figure  (2.1C)  illustrates  this  local  sheet  of  shed  vor- 
ticity. It  should  be  noted  that  vorticity  is  considered  to  be  on  the  helical  surface 
shown    in  Figure  (2.  IB),  the  vertical  displacements  from  that  surface  shown  in  Fig- 
ure (2.1C)  represent  the  strength  of  the  vorticity  at  a  particular  azimuthal  and  radial 
position  [4], 
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Figure  (2.1)  Locwy's  Schematic  Elements  of  Unsteady  Rotor  Plow  Field 
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The  variation  in  this  vertical  displacement  (hence  vortex  strength)  around  the 
azimuth  corresponds  to  the  history  of  the  motion  of  a  given  blade  element  at  a 
fixed  radius.  Variation  of  shed  vortex  strength  in  the  radial  direction  at  any  fixed 
azimuth  angle  is  a  function  of  the  variation  with  blade  span  of  (1)  blade  pianform; 
(2)   amplitude  of  oscillation  of  effective  angle  of  attack,  and  (3)  relative  air  veloc- 
ity.   Since  (from  the  Helmholtz  theorem  )  vorticity  cannot  begin  or  end  in  space,  a 
variation  of  shed  vorticity  in  the  radial  direction  implies  that  there  are  trailing  vor- 
tices at  constant  radii  similar  to  and  inboard  of  the  tip  vortex.  These  trailing  vorti- 
ces have  been  induced  in  Figure  (2.  ID). 
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The  schematic  drawings.  Figures  (2.1  A)  through  Figure  (2.  ID),  Ref.  [4],  indicate 
the  difficulty  of  attempting  to  obtain  a  complete  representation  of  unsteady  rotor 
aerodynamics.  One  means  of  simplifying  the  problem  is  to  consider  extreme  "pitch" 
values  for  the  helix.  For  the  case  of  very  high  inflow  velocity,  u,  in  relation  to  Q, 
and  the  opposite,  where  u  is  very  low  compared  to  Q  This  is  shown  in  Figure 
(2.2). 

When  the   vertical  spacing   between  adjacent  helicai  surfaces  of  shed  vorticity  is 
very  large,  then  we  expect  that  all  shed  vorticity  beyond  a  small  fraction  of  a  rev- 
olution   would  be  too  far  below  the  blade  in  question  to  have  a  significant    effect. 
Under     these  conditions,  it  would  be  sufficient  to  account  for  only  the  attached 
vortex  sheet  within  that  fraction  of  a  revolution,  as  in  Figure  (2. 2A). 

On  the  other  hand,  when  the  vertical  spacing  between  the  adjacent  helical 
surfaces  of  shed  vorticity  is  very  small,  all  the  sheets  of  shed  vorticity  tend  to  pile 
up  on  each  other,  and  the  effect  of  that  vorticity  close  to  the  blade  in  question 
(  shed  by  blade  passes  )  is  of  more  importance  than  that  which  exists  beyond 
a  reference  azimuth  angle  on  either  side  of  the  blade.  This  situation  is  depicted  in 
Figure  (2.2B),  Ref.  [4], 
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lignre  (2.2)  Schematic   Representation  of  Unsteady  Rotor  Flow   Field. 
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With  low  rotor  thrust  values  (  low  collective  pitch   )  the  wake  remains  near 
the  rotor  disk  and  therefore  passes  close  to  the  following  blades.    Thus,  the  shed 
vorticity  sheets  below  the  rotor  disk  must  be  accounted  for  to  correctly  estimate  the 
unsteady  loads.  For  high  inflow  or  high-speed    forward  flight  the  rotor  wake  will 
be  convected  away   from  the  blades.  The  wake  far  from  the  blade  section  will  have 
little  influence,  so    emphasis  is  on  modeling  the  wake  near  the  blade.    In  the  case 
of  low  inflow  this    consists  of  vortex  sheets  that  are  nearly    pianar  surfaces  paral- 
lel to  the  disk  plane  [4] 

Based  on  these  considerations,  a  two-dimensional  model  for  the  unsteady  aero- 
dynamics of  the  rotor  can  be  constructed.      It  is  assumed  that  the  chord,  amplitude 
of  the  oscillation  in  effective  angle  of  attack,  and  relative  airspeed  vary  slowly 
enough  with  span  so  that  what  occurs  aerodynamically  at  one  blade  radius  station 
is  essentially    duplicated  on  either  side  of  it. 

The  assumption  means  that  the  flow  problem  at  a  given  blade  section  is  two- 
dimensional.  Furthermore,  consistent  with  the  idea  that  only  the  vorticity  near  the 
blade  section  has  an  important  effect,  one  may  allow  the  planar  rows  of  vorticity  to 
extend  to  infinity  in  the  horizontal  direction  in  order  to  achieve  mathematical  sim- 
plification. It  should  be  stated  that  this  is  a  first-order  theory,  so  that  such  effects 
as  that  of  the  wake  upon  itself  have  been  neglected.  The  oscillatory  flow  induced 
at  the  airfoil  in  the  free-steam  direction,  for  example,  will  affect  the  airfoil  lift  and 
moment  as  well  as  the  vorticity  shed  at  the  trailing  edge. 
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Since  the  shed  vorticity  is  being  accounted  for  as  perturbations  in  the  main 
How.  this  component  of  induced  velocity  will  be  small  compared  to  the  tree-steam 
velocity  (except  when  the  spacing  approaches  zero)  and  must  react  with  small,  first- 
order  airfoil  displacements.  Because  of  the  physical  unrealities  associated  with  the 
standard  mathematical  model  representation  of  the  problem  (such  as  zero  thickness 
of  wake  and  airfoil,  zero  displacement  of  the  wake  from  horizontal  planes,  and 
infinite  velocities  at  the  core  of  a  vortex),  the  accuracy    of  results  as  the  wake 
spacing  goes  to  zero  is  questionable.  The  effect  of   the  horizontal  component  of 
induced  flow,  therefore,  has  been  neglected  in    Loewy's  analysis. 

Shown  in  Figure  (2.3)  is  the  two-dimensionalized  model  of  unsteady  rotor  flow 
[4].  The  portion  of  the  circular    surface  which  is  determined  by  (1)  a  particular 
blade  radius,  (2)  the  azimuth  angle  on  either  side  of  the  blade  section  (within 
which  the  shed  vorticity  is  of  importance),  and  (3)  the  vertical  distance  spanned  by 
a  given  numbers  of  rows  of  vorticity,    can  be  regarded  as  a  plane,  one  in  which 
the  two-dimensional  unsteady  aerodynamic  problem  may  be  attacked. 

This  is  shown  in  Figure  (2.3  A),  with  an  arrow  indicating  that  the  rows  of 
vorticity  under  the  rotor  disc  extend  to  infinity  Consistent  with  the  idea  that  only 
the  vorticity  near  the  blade  section  has  an  important  effect,  one  may  allow  the  pla- 
nar rows  of  vorticity  to  extend  to  infinity  in  the  horizontal  direction  in  order  to 
achieve  mathematical  simplification.  These  final  modifications  are  shown  in  Figure 
(2.3  B),  which  is  a  two-dimensional  model,  Ref.  [4],  of  unsteady  rotor  aerodynamics 
for  a  single-bladed  rotor  operating  at  low  inflow. 
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Figure  (2.3)  Two-dimensional  Models  of  Unsteady  Rotor  Flow. 
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The  general  approach  used  by  Loewy  is  the  same  as  that  of  Schwarz  [30]. 
Only  the  effects  resulting  from  oscillatory  motion  of  a  thin  airfoil  are  considered, 
and    the  airfoil  itself   is  represented  by    continuous    vorticity    distributed  along  a 
straight    line  [9].    It  may  be  also  noted  that  the  effect  of  viscosity    in  damping  the 
shed  vorticity    could  be  accounted  for  by  considering  only  a  finite  number  of  rows 
beneath  the  reference  airfoil,  or  by  multiplying  the  vortex  strength  by  an  experi- 
mental decay    function  related  to  the  time  through  the  rotational  speed  and  the 
number  of  rotational  cycles  corresponding  to  a  particular  vortex  row. 
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This  has    been  omitted  from  the  analysis  since  vortex  decay    time  is  long 
compared  to  the  rotational  periods  of  most  rotors.  An  extension  to  Loewy  's  work  is 
that  by    Shipman  and  Wood    [29]    in    1970.    The  authors    present  a  method  for 
determining  rotor  blade  flutter  in  forward  flight.  Incompressible  unsteady  aerody- 
namic theory    was  applied  where  shed    vorticity    both  in  and    below    the  plane  of 
the  rotor  is  accounted    for.  This  was  made  possible  by  noting  that  the  damping  for 
rotor  blade  aeroelastic  instability  varies  with  velocity.    Since  the  tangential  velocity 
of  any  blade  segment    varies    with    azimuth,  so  will  flutter  damping  . 

It  is  assumed  at  the  onset  of  flutter  that  oscillations  will  begin    to  build  up 
prior  to  the  blade  reaching  a  critical  azimuth  position,  then  decay    as  the    blade 
moves  beyond  this  point.  This  build  up  and  decay  means  that  the  vorticity  shed 
due  to  the  oscillations  will  be  contained  within  a  double  azimuth  region  on  either 
side  of  the  critical  azimuth  position.  Assuming  this  region  to  be  small  allows  the 
wake  system    to  be  two-dimensional. 

The  lift  deficiency  function  resulting  from  the  two-dimensional  wake  model  is 
compared  with  earlier  results  obtained  for  a  helicopter  in  hover  and  fixed  wings.  In 
the  limiting  case  when  the  advance  ratio  is  very  small,  Shipman  and  Wood's  lift 
deficiency    reduces  to  Loewy's  lift  deficiency  function.  Also  in  the  limit  as  the 
wake  spacing  is  very  large,  this  function  reduces  to  Theodorsen's  lift  deficiency 
function  . 


1  Shipman.  K.  W    ,    and  E.  R.  Wood  "  A  Two-Dimensional  Theory-  for  Rotor    Blade  Flutter 
in  Forward  Flight."     Journal  ol'  Aircraft,  Vol.  No.  8,  Number  12.  pp.    1008-1015. 
December   1971. 
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The  theory  was  applied  to  bending  torsion  flutter  tor  the  tip  segment  of  a  rotor 
blade.  The  results  followed  the  normal  trends  of  having  the  flutter  velocity  decrease 
as  the  blade  center  of   gravitv    was  moved  aft  and  of  having  the  flutter  velocity 
increase  as  the  blade  stiffness  increased  . 

The  previously  shed  wakes  are  destabilizing  so  that  they  reduce  the  flutter 
velocity.  The  buildup  and  decay  produces  wakes  that  are  essentially  centers  of  vor- 
ticity.  Thus  their  position  with  respect  the  reference  blade  will  be  sensitive  to  the 
advance  ratio  which  determines  the  vertical  spacing.  It  was  found    that  as  the 
wakes  were  brought  closer  by  reducing  the  advance  ratio  and/or  inflow  ratio,  the 
flutter  speed    is    decreased  [29]. 

The  phenomenon  of  airfoil  flutter  is  caused  by  energy  extraction  out  of  the 
wind  stream  by  means  of  airfoils  which  are  free  to  oscillate  in  pitch  and  plunge. 
As  described  by  Duncan  [19],  this  effect  can  be  easily  demonstrated  in  a  wind  tun- 
nel if  a  device  is  built  which  allows  to  vary    the  phase  relationship  between  pitch 
and  plunge.  Phase  angles  near  90  degrees  produce  maximum  energy  extraction  so 
that  the  device  operates  as  a  flutter  engine.  McKinney  and  DeLaurier  [20]  demon- 
strated that  this  effect  can  be  used  to  convert  wind  energy  into  mechanical  energy. 
Scherer,  McKinney  and  Delaurier  used  Theodorsen's  flat-plate  theory  for  perfor- 
mance estimates  of  the  flapping  foil  propulsor  and  of  the  wingmill,  respectively. 

It  was  shown  by  W.  Schmidt,  that  optimum  energy  extraction  occurs  for  har- 
monically flapping  airfoils  upstream  of  a  stationary  airfoil.  There  is  a  need  for 
aerodynamic  analysis  methods  which  enable  systematic  studies  of  the  effect  of  air- 
foil geometry,  incidence  angle  variation,  amplitude  of  oscillation  and  airfoil  interac- 
tion on  the  aerodynamic  forces  due  to  incompressible  flow  over  airfoils  or  airfoil 
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combinations  executing  pitch  or  plunge  motions.  A  panel  method  was  developed  in 
order  to  do  parametric  studies  of  these  effects.    In  the  past,  a  number  of  investiga- 
tors have  solved    the  steady    flow    problem  using    source  and    vortex  paneling,  the 
most  prominent  ones  being  Hess  and  Smith  [21].  A  few  authors  have  extended  this 
approach  to  the  case  of  unsteady  motion  of  single  airfoils,  notably  Basu  and  Han- 
cock [22]  and  Kim  and  Mook  [23]. 

Basu  and  Hancock  [22]  presented  a  numerical  method    for  the  calculation  of 
the  pressure  distribution,  force,  and  moments  on  a  two-dimensional  airfoil  undergo- 
ing an  arbitrary  unsteady  motion  in  an  inviscid  incompressible  flow.  The  method 
was  applied  to  (i)  a  sudden  change  in  airfoil  incidence,  (ii)  an  airfoil  oscillating  in 
high  frequency  and  (iii)  an  airfoil  passing  through  a  sharp-edged  gust. 

Teng  [24],  developed  a  computer  code  for  the  numerical  solution  of  unsteady 
inviscid  incompressible  flow  over  an  airfoil.  Pang  [25]  extended  the  work  of  Teng 
and  developed  a  computer  code  to  study  unsteady  airfoil  interference  effects.  The 
technique  uses  the  well  known  panel  methods  for  steady  flow  and  extends  it  to 
unsteady  flow  by  accounting  for  the  continuous  shedding  of  vortices  into  the  trailing 
wake. 
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III.  DRAG  REDUCTION  DUE  TO 

BLADE    FLAPPING 

(  KATZMAYR  EFFECT  ) 

A.   INTRODUCTION 

In  this  chapter,  a  discussion  will  first  be  presented  of  the  physics  of    Happing 
airfoils  which  explains  the  bird's  ability  to  generate  a  forward    thrust    by    means  of 
wing  flapping.    This    will    be  followed    by    a  description  of  two   methods    capable  of 
predicting  the  pressure  distributions,  lift,  pitching  moment,  and  drag  or  thrust  on  airfoils 
which  execute  harmonic  oscillations  in  plunge  or  pitch. 

The  first  method  is  based  on  the  generalization  of  the  well-known  steady-state 
panel  method  to  the  case  of  unsteady  airfoil  motion.  It  is  based  on  the  assumption 
of  inviscid.  incompressible  flow.  The  second  method  is  based  on  the  classical  The- 
odorsen  theory  for  oscillating  fiat-plate  airfoils  in  incompressible  inviscid  flow.  The 
unsteady  panel  method  permits  the  investigation  of  finite-thickness  airfoils  of  arbitrary 
geometry  which  may  execute  rather  general  unsteady  motions  (  ramp  and  oscillatory  ). 
In  contrast,  Theodorsen's  analysis  is  restricted  to  thin  airfoils  ( flat-plates )  and  to  small 
amplitude  oscillations.  The  chapter  ends  with  a  comparison  of  results  computed  by  both 
methods. 
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B.    THE  PHYSICS  OF  FLAPPING  AIRFOILS   IN   LOW-SPEED   FLIGHT 

The  phenomenon  of  thrust  generation  due  to  wing  flapping  or  due  to  steady  flight 
in  a  sinusoidal  gust  can  be   explained  in  a  rather  elementary  way.     Consider  the  latter 
case  of  flight  through  a    sinusoidal  gust,  as  shown  in    Figure  (3.1).  The  relative  wind 
vector  changes  direction  between  a  maximum  positive  and  negative  angle  of  attack.  The 
resultant  aerodynamic  lift  is  perpendicular  to  the  instantaneous    wind  vector.    If  we  as- 
sume the  viscous  drag  to  be  small,  then  it  is  easilv  seen  that  the  sinusoidal  gust  gen- 
erates a  sinusoidally  varying  thrust. 

A  similar  effect  is    generated    if  the  airfoil  executes  a  sinusoidally  varying  plunge 
motion  (  wing  flapping )  about  an  otherwise  steady  flight  conditions.  This  explanation 
of  thrust  generation  due  to  wing  flapping  or  due  to  flight  through  a  sinusoidal  gust  was 
first    advanced    by    Knoller  [10]  and  Betz  [11].    Its  experimental  demonstration  was 
accomplished  by  Katzmayr  [7]  and    is   therefore  generally    referred   to   as   Katzmayr 
effect  l . 

Since  it  is  well  known  that  every  change  in  airfoil  incidence  causes  the  shedding 
of  a  starting  vortex,  a  more  precise  explanation  of  the  Katzmayr  effect  has  to  take  into 
account  the  continuous  vortex  shedding  due  to  wing  flapping  or  due  to  flight  through 
a  sinusoidal  gust. 


1  Platzer.  M.  F..  Neace.  K.  S  ,  and  Pang.  C.  K.. "Aerodynamic  Analysis  of  Flapping  Wing    Propulsion", 
page  1.  AIAA  paper  No.  93-0484.  31st    Aerospace  Sciences  meeting.  Reno.  NV.  January    11-14.   1993. 
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Figure  (3.1)     Propulsive  force  on   plunging  airfoil 
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Similarly,  the    quantitative    prediction    of  the   induced    oscillatory    forces  (  lift  and 
thrust )  had    to    await    the  development  of  a  complete  unsteady  aerodynamic  theory. 
Such  a  theory    was  first  accomplished  in  the  1930'  s  by  Theodorsen  [3]  in  the  United 
States  and  by   Klissner  [16]  in  Germany.    It  can  be  found  in  the  textbooks  on  aeroelas- 
ticity,    for   example  the  book   by    Bisplinghoff,  Ashley,  Halfman  [41].    Therefore  only 
the  final  formulas  will  be  given    here    and    the  reader    is    referred    to    those  texts  for 
additional  details.    More    recently,    a    second    approach    to  the  determination    of   the 
unsteady  aerodynamic  forces  on  oscillating  airfoils  became  feasible  with  the  advent  of 
high  speed  computers  and  the  development  of  efficient  numerical  solution  methods.   The 
fundamental  building  blocks  of  this  method  are  presented  in  the  next  two  sections. 

C.  PANEL  CODE  FOR  STEADY  INCOMPRESSIBLE  INVISCID  FLOW 

In  potential  flow  theory,  the  flow  field  around    an  airfoil  may  be  represented  by 
the  velocity  potential.    Considering  contributions  from  the  free  stream  flow  and  the 
source  and  vorticity  distribution,  the  total  potential  may  be  constructed: 

(D    =    (Dp    +    <DS    +    <DV  (31) 

where 

<Dp  =     U  [  x  cos  (  a  )  +  y  sin  (  a  )  ] 


q(s) 

cd  =-f^eds 


<D    =   p^Iin(r)ds  (3.2) 

T(s) 


m 
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The  source  distributions  q(s)  vary  from  panei  to  panel,  while  the  vorticity  strength 
y  (s)  is  assumed  constant  for  ail  panels.    The  value  of  representing  the  flow  past  an 
airfoil  by  surface  singularity  distributions  lies  in  the  fact    that  these  singularity    distri- 
butions automatically  satisfy  Laplace's  equation,    the  governing  flow  equation  for  invis- 
cid  incompressible  flow: 

0>     +0      =0  (3.3) 

\  \  v  y 

Since  Laplace's  equation  is  a  linear  homogeneous  second  order  partial  differential 
equation,  the  superposition  principle  used  in  Equation  (3.1)  holds.  The  boundary  con- 
ditions include  flow  tangency  at  control  points  (mid-points  of  panels)  and  the  Kutta 
condition  at  the  trailing  edge,  requiring  equal  tangential  velocities  for  the  first  and  last 
panels.  By  evaluating  the  integrals  along  the  airfoil  surface,  the  potential  may  be  de- 
termined at  any  point  in  the  flow  field.  Each  point  is  defined  at  a  radius  (r)  and  angle 
(6)  from  a  chosen  reference  point  on  the  airfoil. 

The  reference  point  in  this  study  is  the  leading  edge.  The  airfoil  is  represented 
by  a  number  of  defined  points,  called  nodes  [34].  More  points  produce  greater  resolution 
and  better  accuracy      One  hundred    to  two  hundred  points  are  usually  sufficient,  with 
the  larger    number   used  for  more  complicated  airfoil  shapes  or  more  involved  calcula- 
tions  The  lines   connecting  these  nodes  are  the  panels.  There  are  n+1  nodes  with  first 
and  last  node  overlapping. 
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■igure  (3  2)  depicts  the  panel  geometry.  Numbering  starts  at  the  trailing  edge,  then 
resses  along  the  lower  surface,  leading  edge,  and  upper  surface,  and  ends  at  the 
ng  edge.  The  unit  normal  vectors  (n,)  are  perpendicular  to  the  panels  and  directed 
aid   from  the  airfoil   surface    The  unit  tangential   vectors  (t,)    are  parallel   to  the 
Is  and    the  positive  direction  is  defined     with   increasitm  numbering  (  n  (o  tiH  ) 
panels  may    vary    in  length,  with  the  exception    of  the  llrst  and    last  panels,  which 
be  equal   in    order  to  use  the  Kutta  condition    at  the  trailing  edge. 


^      I 


n-i    n 


n+L 


Figure  3.2  Panel  Method  Geometry 
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The  required    input    consists  ot    the  number  of    nodes  on    the  airfoil  surface,  the 
coordinates  of  the  nodes  referenced  from  the  leading  edge,  and  the  angle  of  attack  in 
degrees.  The  program  produces  normalized  velocities  and  pressure  coefficients  at  each 
control  point  as    output.    The  use  of  influence  coefficients  leads  to  a    straightforward 
procedure  for  programming  the  equations. 

An  aerodynamic  influence  coefficient  is  defined  as  the  velocity  induced  at  a  field 
point  by  a  unit  strength  singularity  distribution  on  one  panel.  For  the  two-dimensional 
steady  flow  problem,  the  following  influence  coefficients  are  needed  [341 

Ann  .     normal    velocity  component  induced  at  the  i  ;1  panel  control  point  by  unit 
strength    source  distribution  on  the  il     panel 

1  r       l 

An  =   —    sin  (8.-9.)  In— —  +  cos(0.  -0.)B-     -        i*j  (3.4) 

i]        ~>K  i       r         r  x        '      lJ 

»J 

=  1/2  i=j 

Aljj  :  tangential  velocity  component  induced  at  the  i  panel  control  point  by  unit 
strength  source  distribution  on  the  i      panel 

A1.  =    —    sin(0.-e.)p..-cos(0-G)ln-^      ,      i  *  j  (3.5) 

ij        2%  '        '      "  '        '  r 

ij    _ 

=0  ,       i=j 

B",,  .     normal    velocity  component  induced  at  the  i      panei  control  point  by  unit 
strength  vorticity  distribution  on  the  j      panel 

1                               r  ■     i 
Bn.  =  cos(6.-0.)ln-         -sin(0.-0.)p..      ,       i^j     (3.6) 

ii        ",ti  ''  r  'in 

u 

=  0  .      i=j 
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Bln       tangential    velocity  component  induced  at  the  i      panel  control   point  by  unit 
strength  vorticity  distribution  on  the  j      panel 

Bl    =  cos  (0. -0.)  B.. -i- sin  (9.-9.)  In  —        .  i  =j  (3  7) 

ij        2k  '        '      !|  '        i  r 

ij    - 

=   1/2  .1=1 


where  the  geometrical  quantities,  depicted  in  Figure  (3.3)  are  defined  by 

x .  -t-  x 
li-i 
r     =      (xm   -  x.)     -  (vm  -  v  )     .        \m    = 

11  !]  J         1         "  J  1  "> 

V  ■  -  V       .  V       i  -  v 

vm    =  0.  -  arc  tan  (j.8) 

"      '2  '  x .     ,  -  x .  j 

i  +  ] 


p..  =  arctan  I  I  -  arctan 

lJ  Um  -  x.     ,  J  I  xm.  -  x.  J 

v       i       j  +  l y  l       J 

The  first  boundary  condition  requires  flow  tangencv  at  the  control  points 

(Vn  ),  =  0  .        i  =  1,  2, n  (3.9) 

In  terms  of  influence  coefficients  (with   V     =   1), 


V     Anq.    +y  V  Bn^sin(ct-e  )    =  0  ,      i=l,2 n         (3.10) 

j  =  1  j  =  1 

The  second  boundary  condition  is  the  Kutta  condition,  which  states  that  the  pressures 
on  the  lower  and  upper  panels  at  the  trailing  edge  must  be  equal  if  the  flow  is  to  leave 
the  trailing  edge  smoothly. 
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Figure  (3.3)    Relationship  Between  Geometrical  Quantities 
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Using  Bernoulli's  equation,  the  pressure  coefficient  can  be  expressed  as 

r  loin 

C     =  I  -  -  (3.1 

P  I         2  \     V       J 

pV  ^  • 

2 
The  pressure   equilibrium  also  implies  equal  velocities  for  incompressible  How    Since 

the  normal  velocities  are  taken  to  be  zero,  the  boundary  condition  may  now  be  stated  as 

(  V  ),     =  -    (  V1  )„  (3.12) 

wheic  the  negative  sign  is  strictly  due  to  the  adopted  convention  of  positive  tangential 
velocities  in  the  direction  of  increasing  node  numbering.  Since  the  How  is  positive  to 
the  light,  the  panel  downstream  of  the  front  stagnation  point  will  have  negative  values 
for  computational  purposes  only  It  is  important  to  note  that  not  all  the  lower  surface 
panels  have  a  reversed  sign,  only  those  downstream  from  the  stagnation  point.  This  is 
especially  significant  for  non-  symmetrical  airfoils  or  any  airfoil  at  an  angle  of  attack. 
In  terms  of  influence  coefficients,  the  normalized  equation  becomes 


it  n  n  n 


.  t 


V     A..q.    -y  V  B..  -cos  (a -8.)   =    Y     A   q.     -  y  y  B..  -  cos  (a  -  0,)      (3.13) 

t—t    I       ll  'l  I        '   Z-i       U  I  Z-i  il  M  '   L-i       II  \' 


t 


1 


i.l  M  i         *-*      u 

i  =  i  i  =  i  1  =  1  i  =  i 


Equations  3.10  and  3  13  represent  a  linear  algebraic  system  of  (n+1  )  equation  and 
(nH)  unknowns.    The  unknowns  are  the  source  strengths    which    vary    from  panel    to 
panel  (    qt  qn  )  and  the  vorticity  strength    y. 
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Expanding  and  rearranging  Equation  5.10  for  an  example  airfoil  of  n  =  73  noaes  and 
)aneis  results  in 


A!lql  "  A?2(h  +  ■■■ 

-•/(b;',  -b?.  - ... 

-  B'? 

_,;    =    -Mina 

-e.) 

A?,q,  ~  A_q:  - 

rydJ!!,  -Bt-... 

-b;1 

_,,    =    --in  (a 

-e,) 

\  n      ■                 v  :! 

...  -y(bL\     -  b_, 

- 

..  -B2,  _,;    - 

-  -in 

(3.14) 


The  equations  now  readily  lend  themselves  to  soiution  in  matrix  form  Recasting 
with  simpler  notation,  the  A;\,  terms  (coefficients  ot  q,  )  mav  be  renamed  a„  and  the 
sum  of  ail  Bn  ,,  terms  in  parentheses  (coefficients  of  y  )  renamea    a,  nr,;  .  where  i  =   1. 

2.  n    and    i  =  1  .  2,  n     The  terms  on  the  right  side  of   the  equation    may 

be  renamed  b,  .     The  (n+1)  equation,  or  in    this  example,  the  74lh  equation    comes 
from  Equation    3.13  in  a  similar  manner: 

(A^  !  -AL,   l)q,  +  (A1,  ;  rA;3i2)q,...  -  (A1.   _,  -  A_)q_, 

-y[  (  Bl, A-  Bl-,j  )  -    (  Bl12-  Bl -,.:)  +  -  (   BlL73+  Bl -, -3  )  ] 

=    cos  (a -9])-    cos(a-0-73)  (3.15) 

The  coefficient  q,  may  be  renamed  a-4  ,  All  of  the  Bl  terms  in  the  brackets  together 
form  the  coefficient  of  y,  now  renamed  a  74  74  .  The  entire  right  side  of  the  equation 
constitutes  the  new  term  b-74 
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Finally  expressing  this  system  can  be  written  in    concise  matrix  form 

i 
al,  I      a!.2 Vri+l       |q,  bl 

a:.i     a:.: a2.n+i     |q2  b: 

(3.16) 

a     .      a     , a  a  b 

n.  1         i.  J  n.  n  +  I        Mn  n 

a  a     a  ,   i  v  b 

n+  I,  1  n  +  I,  n  +  1    L  i  _  n  +  I 

This  system  is  solved  in  the  program  using  a  Gaussian  Elimination  subroutine.    With 

the  values  of  the  u,    and  y  known,  the  velocity  at  each  panei  control   point  may  be 

calculated 


V    =    V     A  q     +cos(a-0)     ,    i  =  1,  2,  n.  (3.17) 

i  Lu         ii  l\  v 


n 

t 

il  'I 
J  =  I 

The  total  velocity  is  equal  to  the  tangential  velocity  due  to  taking  the  normal  velocity 
to  be  zero. 
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D.    UNSTEADY  PANEL  CODE 

The  extension  of  this  method  to  an  airfoil  experiencing  :.  rapid  change  in 
.ingie  oi  attack  can  be  accomplished  as  follows,  it  is  important  to  recaii  Prandtl's 
eariv  flow  \  isuaiization  experiments  which  demonstrated  the  generation  of   a 
starting  vortex  off  the  trailing-edge  whenever  the  angle  of  attack  was    changed. 
Hence  continuous    changes  in  angle  oi  attack  oroduce  a  continuous  shedding  of 
vorticitv  into  the  trailing  wake    This  shedding  can  be  expiaineu  by  the  Heimholtz 
vortex  theorem  which  requires  that  any  change  in  the  circulation  around  the  airfoil 
must  be  matched    by  the  appearance  ot  an  equal  counter  vortex  at  the  trailing 
edge  in  order  to  achieve  constancy  ot  the  total  circulation, in  me  flow  field. 

The  vortices  shed  from  the  trailing  edge  move  with  the  fluid  particles  of  the 
surrounding  fluid  and  hence  are  swept  downstream  with  a  speed  essentially  equal 
to  the  free-stream  speed.  Therefore  these  vortices  will  stav  close  enough  to  the 
airfoil  for  a  finite  time  to  influence  its  pressure  distribution.  Any  unsteady  airfoil 
theory  therefore  must  describe  this  vortex  shedding  process.  It  is  this  feature 
which  distinguishes  unsteady  airfoil  theory  from  its  steady  counterpart'  The 
above  considerations  suggest  that  the  airfoil  can  be  modeled  by  similar  source 
and    vortex    distributions  as  in  the  steady    case. 


Platzer.  M.  F.   Class  Lecture  Notes.  Naval  Postgraduate  School.  Monterey 


California.   Sen.    i 993 
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Hence    source  panels  are  again  placed  on  the  airfoil  surface  together  with  a 
single  vortex  distribution  which  is  the  same  for  each  panel.  An  arbitrary  change  in 
angle  of  attack,  then  is  subdivided  into  small   step  changes  such  that  a  starting 
vortex  is  shed  into  the  wake  each  time  a  step  change  occurs  and  the  flow  tangencv 
and  Kutta  conditions  are    enforced  at  each  time  step. 

In  addition  to  the  n  airfoil  panels  awake  panel  is  assumed  to  be  attached  to 
the  trailing  edge  such  that,  after  each  time  step,  the  vorticity  ot  the  wake  panel  is 
assumed    to  be  concentrated  into  a  single  point    vortex  which    detaches  from  the 
airfoil  and  starts  moving  downstream  with  the  free-stream    velocity.     \s    one  wake 
panel  detaches  from  the  airfoil   a  new  one  is  created  and  a  wake  of  point  vortices 
is  created.  As  in  the  steady  problem  n  unknown  source  strengths  and  one  unknown 
vorticity  strength  on  the  airfoil  are  introduced. 

In  addition,  there  is  now  the  unknown  vorticity  strength   of  the  wake  panel 
whose  length  and  orientation  are  also  unknown.  Hence  there  are  three  additional 
unknowns    which    need  to    be  determined    by    a  proper  set  of  equations;    (a)  the 
flow  tangencv  conditions  provides  n   equations,  (b)  the  Kutta  condition  provides  an 
additional  equation,  (c)  the  Helmholtz  theorem  provides  an  equation  for    the  vortic- 
ity strength  on  the  wake  panel. 

Therefore  two  more  equations  are  needed  to  determine  the  length  and  orienta- 
tion of  the  wake  panel.  These  are  obtained  by  making  the  following  assumptions,  i.e.. 
(d)  the  wake  panel  is  oriented  in  the  direction  of  the  local  resultant  velocity    at  its 
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nid-point.  ;e)  the  length  of  the  wake  panel  is  proportional  to  the  magnitude  01'  the 

local  resuitant  velocity  at  its  mid-point  and  to  the  nme-siep      The  x.una  condition 

and    the  last  two  conditions  are  nonlinear    equations    and    therefore  necessitate  an 

iterative  solution  procedure     Figure  (3  4).  from  Ref    '251.  shows  the  panel  metiiod 

representation  for  unsteadv  flow 
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Figure   (3.4)    Pnncl    Methods    Representation    for   Unsteady    Flow 
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E.    HIEODORSEVS  ANALYSIS 

As  stated  before.  I'heoaorsen's  theory  for  osciilating  flat-piate  airtbiis  in  inviscid 
incompressible  flow  is  well  documented.  Therefore.  e  summarize  here  :miy  the 
final  formulas.  The  iift  on  an  airfoil  oscillating  at  frequency  co  is  given  by 

L  ==  -ipo:  b';  :   Lh  .  (h,y'b)     -  [  L,x    -  (    1.2        i  )    U     i    a   J  i  3.18) 

\  here 

Lh    =   !   -  i  2i    COO  )      k  . 

La    =    1/2    -  ( i  /  k  ).  [  1  +  2  C(k)  ]  -     (  2  /  k2  )  C  (k). 
k  =    co  b  /  U       (    b  is  the  semi-chord    and  "a"  is  the  elastic  axis  ) 
C(k)    is  the  lift  deficiency    function    and  is  mven  bv 


H,,2)(k) 

C(k)   =  =     F(k)    -  i  G(k)  (3.19) 

H  -l  (k)  -.HU)(k) 


where    H"n(k)    is  the  Hankel  function  ot  second  kind    and  is  given  in    terms 
ot    Bessei  functions  as  :        H~n(k)     =    Jn(k)    -    Vn(k) 
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F.    GARRICKS  PROPULSIVE   FORCE  ANALYSIS 

Garrick  [5]  applied  Theodorsen's  theorv  to  determine  the  thrust  and  drag 
generated  by  an  airfoil  which  executes  a  sinusoidal  motion  in  pitch  or  in  plunge. 
He  found  the  following  formula  for  the  average  propulsive  force  : 

Px  =  TtpbeT    a,  h~  -  (a.  -bja"  -2(a    -b)a  (h  (3.20) 

where         aj    =  F~  -  G" 

a2         b:  !  a,  [  1  k~  ■*-  (  1/2  -a)2  J  +1/4  -  ( 1/2  -  a  )  F  -  G  /  k  I 
a4  =  b   J   aj  [  -  sin  (  (p2  -  (0o  )  +    (1/2  -  a  )  cos  ((p2  -■  (0O)  ] 
-  (F/2)  cos  (  (p2  -<p0)  +  (G/2) sin  ( (P:  -  (Po)  ! 
b2   =    b2[  -  a/2    -  F/k2   +  (1/2  -  a  )  (G/k)  ] 
b4     =  (b/2)  [  (1/2  +  G/k )  cos  (  (p2-  ip0  )  +  (F/k)  sin  (  (p2  -  cp0)] 
For  pure  plunge,  the  average  propulsive  force  is  reduced  to 

Px  =  rcpboV   F2(k)  +G2(k)  ]  (3.21) 

For  pure  pitch,  the  expression    for  the  average  propulsive  force  is  given  by  : 

3     "*     2 
?a  =  7tpb  (o~a~  [a^ -b.]  (3  22) 

where  a2    and    b2  are  defined  in  Equation  (3.20). 
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The  expression  for  the  average  drag  in  Equation  (3.20)    of  Garrick's  reDort 
"o?  1 5], was  found  to  be  m  error    The  correct  cxoression  as  given  in  Garnck's 
subsequent    review    [32]  is 

i     -     I        ■ '  1                        i              1                  F        '  1             G  i 
D  =  A<C  -  v  -F  -X  -v  ,3  23) 

vv  here 

C"  =  F~    -  G"     and     A  =  -  z .  p.  co~.  at>~  b"" 

G.  COMPARISON   OF  RESULTS 

in  a  recent  thesis.  Riester  [  43  ]  presented  a  uetaiied    comparison    with 
Theodorsen's  predictions  ot  the  lift  and  pitching  moment  comDuted  with  the  un- 
steady panel  code.    Generally,  good  agreement  was  found,  thus  lending  credence 
to  the  reliability    of  the  panel  code.    In   the   present    work,  we  are  interested  in 
the  prediction  of  the  thrust  and  drag  forces  generated  by   an  oscillating  airfoil. 
Calculations  were  performed  to  ascertain    the  panel  code's  capability      Figures 
(3.5)  to   (3.7)  show  the  comparison   between   Garnck's  predictions  and  the  oanei 
code  results.    These  figures  present  plots  of   the  average  propulsive  force  coef- 
ficient vs  the   non-dimensional  plunging  amplitude  of  oscillation.  HBAR  (FffiAR 
=  h/'c  ),  for  three  reduced  frequencies,  k  =  0.1  ,   I  and  2.    The  two  methods  are 
in  good  agreement  for  small  plunge  amplitudes.     As    might  be    expected,  the 
agreement  starts  to  deteriorate  with  increasing  plunge  amplitude.   The  panei  code 
was  applied  to  an    \ACA  0009  airfoil  using   100  panels. 
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Figure  (3.5)  Comparison  of    Panel    Code  Results  with  Garri 
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H.   ESTIMATED    POWER  REDUCTION    Dl  E   TO   HHC  FOR  OH-6A 
HELICOPTER 

The  ultimate  objective  of  this  research  was  to  determine  the  effect  or  unsteady 
aerodynamics  on  the  measured  power  values  for  the  OH-6A  helicopter  We  are 
especially  interested  in  knowing  whether  these  effects  couid  produce  the  power 
changes  recorded  during  the  NASA  -  Army  open  loop  flight  tests  in  the  flight 
regime  from  hover  to    100  knots. 

In    order  to  estimate  the  change  in  torque  due  to  propulsive  force  (  foot- 
pounds per  blade)  due  to    HHC.  the  blade  ^pan  was  divided  into  sections     As- 
suming hover  conditions,  the  reduced  frequency  at  each  blade  section  was  cal- 
culated   The  panel  code  was  used  to    compute    the  propulsive  force  coefficient. 
Cd,  at  each    blade  section  assuming  pure  plunge    motion    with  a  plunging    am- 
plitude of   1".    The  blade  deflections,  obtained  from  Wood   et  al  [27],  were  then 
used  to  estimate  the  real  value  of  the  plunging  amplitudes  at  each  blade  segment. 

Finally,  the  net  blade  torque  could  be  calculated  by  considering  the  propulsive 
force  at  each  blade  section  and  multiplying  this  value  times  the  rotor  radius,  then 
summing  over  the  length  of  the  blade.  Table  3.1,  from  Ref  [42],  shows  the 
OH-6A  helicopter  parameters  that  were  used  to  estimate  the  propulsive  forces. 
The  most  important  parameters  are:  number  of  blades.  0=4.  rotor  radius  (ft), 
R  =  13.2.  blade  chord  (ft),  c  =  2b=  0.57,  rotor  speed  (rad/sec).  Q  =49.2.  and 
helicopter  weight    (lb),  G  =  2550. 
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TABLE  3.1  THE  OH-6A  HELICOPTER  ROTOR  PARAMETERS 

Rotor  radius  (ft)  13  2 

Number  of  blades  4 

Rotor  speed  (rad/sec)  50  26 

Blade    chord  (ft)  0.57 
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Table  \3. 2)  shows,  at  each  blade  section  r/R.  the  reduced  frequency.  k(r). 
local  speed.  UCr),  the  computed  average  propulsive  force.  P.  ,  (r)  ana  the  torque. 
TM1',i)  aue  to  the  propulsive  force  at  each  blade  section,  humming  the  esti- 
mated torque  at  each  blade  section,  the  resultant  torque  due  to  the  pronuisive  force 
per  blade  ^as  found  to  be  8  1 1  ft-lb  assuming  a  pure  piungmg  motion  with 
constant  piunging  amplitude  of  I  "  at  each  blade  section.  Wood  et  ai  [27].  ap- 
plied Garrickfs  equation  and  considered  the  effect  oi  the  second  to  twelfth  har- 
monic loads.  He  found,  for  2P  harmonic  (  ITHC  :n  '  330  degree  phase  at  60 
knots),  that  the  power  gained  is  about  11.3  horsepower  For  more  detaii.  the 
reader  is  referred  to  table  2  oi  Ref.  [271  Figures  3  3  to  -•  !  1  i.'tow  the  blade 
deflections  for  several  harmonics  used  by  this  reference.  hown  in  figure  (3  12) 
is  the  average  drag  coefficient.  Cdav.  \ersus  reaucea  frequency  as  obtained  with 
the  panei  code  for  a  plunging  amplitude  of  1"  The  airfoil  used  was  N'ACA  0015 
(OH-6A  airfoil).  The  propulsive  force.  Pxll ,  is  given  by 

PNh  =  0.5.  p.  (J2    C.dr.  Cd  (3  24) 

and  the  torque.  Txh.    at  each  blade  section  is  given  by: 

where    U  is  the  local  speed    at  each  section  ,    Cj  is    the  propulsive  force  coef- 
ficient as  obtained  by  the  unsteady  panel  code  and    r    is    the  radial  distance  oi 
each  blade  section    For  the  OH-6A.  the  angular  frequency  .  CI  =  50.26  rad/sec, 
therefore  the  power  per  blade  is  found  to  beO.7-1  horsepower  as  given  bv  : 

P  =  Txh    Q  .  (3.26) 
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TABLE  3.2  LOCAL  REDUCED  FREQUENCY.  K?(r).   PROPULSIVE 
FORCE.  Pxh(r)  AND  TORQUE.  Txh  (r) . 


r/R  k(r/R) 


U(r/R) 


Pxh(r/R) 


T,h(r/R) 


[  "  ]   i  [ 


ft  /  sec 


,r  Ibf 


rt-lb 


0. 1 

0.2 
0.3 
0.4 

0.5 

06 

0.7 
0.8 
0.9 


0  863 
I ) 43  1  5 
0.2877 
0.258 
0.173 

0.144 
0.1233 

0.1079 

0  0959 


0.0618 


0  0798 

0.9530 

0.108 

0  118 

0.127 

0.134 

0.139 

0  144 

0.122 

0  26 
0  44 
0  64 
0.86 

1.09 
1.32 
1.56 
1.81 


VT,    -81    tt-lb 
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IABLE  3.3  :  (  d.,v    versus  REDUCED  FREQUENCY    K"    at   hbar  =  I 


kp       2  k 

Cdnv 

0. 1 

-0.0004937 

0.2 

-0.0018445 

0.3 

-0.003588 

0.4 

-0.0056899 

0.5 

-0.008108 

0.6 

-0.01082 

0.  7 

-0.01363 

0.8 

-0.01724 

0.9 

-0.02057 

1 

-0.02457 

1.1 

-0.029702 

1.2 

-0.033156 

1.3 

-0.0379 

1.4 

-0.04291 

1.5 

-0.04822 

1.6 

-0.05382 

1.7 

-0.05972 

1.8 

-0.06592 

1.9 

-0.07241 

2 

-0.079203 

2.1 

-0.08633 

1 
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IV.  ANALYSIS  OF  WAKE 
INTERFERENCE   EFFECTS 

A.  INTRO  DICTION 

The  results  of  Chapter  III    show    that    propulsive  forces,    generated    solely 
by     plunge    oscillations    of   an    airfoil  are  insufficient  to    explain    the  power 
reduction  observed  and    measured    on  the  OH-6A  helicopter     Therefore,  in  ihis 
chapter,  we  examine  whether  the  power  reduction  might  be  caused   by   favorable 
wake  interference  effects  due  to  the  wakes    shed    from    preceding  blades. 

To  this  end.  the  extension  of  the  unsteady  panel  code  to  the  case  of  two 
airfoils  [25]  was  adopted  and  applied  to  the  study  of  the  wake  interference  effects 
between   two  helicopter  blades.  Wake  interference  effects    were  studied   in  the 
1950's    by    Loewy    [4]  by  extending  Theodorsen's  analysis.     Loewy  treats  an 
infinite  number  of  layers  of  shed  vorticity  placed  beneath  the  rotor  at  a  given 
wake  spacing  and  developes  a  modified  lift  deficiencv    fiinction      His  interest 
was  in  the  possibility  of  wake-induced  blade  flutter,  which  restricted  his  analysis 
to  the  determination  of  lift  and  pitching  moment     changes  due  to  wakes  shed 
from  preceding  blades.  A  closer  inspection  of  Loewy's  paper  suggests     that  it 
should  be  possible  to  use  the  same  approach  to  examine  drag  and  propulsive 
force  changes  on  the  reference  blade  due  to  wakes  shed  from  neighbouring  biades. 
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Further,    it    is    desirable  to  carry  this  one  step  further  and  investigate  the 
ease  of  a  single  wake  shed  from  the  preceding  blade  which  interferes  with  the 
wake  shed  from  the  reference  blade      This    case    allows    a  direct  comparison 
with  the  unsteady    two-foil    panel  code. 

Also,  there  are  regimes  of  helicopter  flight  where  oniy  the  first  laver  of 
shed  vorticitv  is  of  consequence.    In  the  following  sections  we   first  describe  the 
two-foil  panel  code  and  Couch's  modification  of  the  Loewy  analysis     This  is 
followed  bv    a  presentation  and  comparison  of    the  major  results. 

B.   UNSTEADY  TWO-FOIL  PANEL  CODE 

The  extension  of  the  single  airfoil  code  (  previously  described  )  to  the  case 
of  two  airfoils  requires  no  new  building  blocks.  However,  the  two-foil  analysis 
requires  the  introduction  of  five  frames  of  reference,  namely  two  moving  local 
frames  of  reference  which  are  attached  to  the  two  airfoils,  two  frozen  local  frames 
of  reference,  and  the  inertial  frame  of  reference.  Furthermore,  it  requires  the 
satisfaction  of  the  two  Kutta  conditions  which  are  coupled  non-lineariv  The 
solution  requires  an  iterative  procedure  to  compute  the  two  vorticitv  distributions. 


A  detailed  development  ol~  this  theory  is  derived  and  presented 
in  a  recent    NFS  thesis  bv   Mark  Couch.  Ref.   1401 
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It  also  necessitates  the  creation  of  a  subroutine  which  transforms  ail  coor- 
dinates m  either  or'  two  local  frames  of  reference  to  the  global  frame  of  reference. 
Finally,  it  requires  the  extension  of  the  influence  coefficient  concept  to  include 
the  effect  of  the  second  airfoil  with  its  own  wake  and  it  requires  the  introduction 
of  an  additional  influence  coefficient,  namely  that  on  die  wake  element  jue  to 
the  wake  element  from  the  other  airfoil.  !f  each  airfoil  is  modelled  with  X 
panels,  then  this  produces  a  2Nx  <  N'+S  )  matrix  which  is  subseauentiy  solved 
by     Gauss  elimination 

Pang  [25]  was  only  able  to  verify    this    code    by    .ompanng  me  computed 
pressure  distributions  with  previous  work  by    Giesing  [531      Hence  it  was  nec- 
essary to  further  evaluate  the  code  by  comparing  its  output  with  otner  known 
solutions   before  applying  it  to  the  oscillatory   blade  interference  problem.  These 
comparisons  are  given  in  section  C. 

The  main  parameters  used  in  this  code  are:   the  airfoil  type,  number  of 
paneis.  N.  the  relative  locations  of  the  two  airfoils,   initial  angle  of  attack,    am- 
plitudes of  oscillation  (pure  pitch,  pure  piunge  or  both),  frequency  of  oscillation 
for  each  airfoil,  the  pivot  points  (for  pitch  cases),  rise  timet  for  ramp  motions  ). 
These  input  parameters  are  in  the  file  "fort.  1" 

This  program  produces  an  extremely  large  amount  of  output  to  the  screen.  It 
is  convenient  to  write  the  output  to  a  file  on  a  L'nix  based  machine  during 
program  runs.  The  output  file  can  get  very  large  for  long  program  runs  (when 
the  time  history    is  long). 
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For  this  reason,  the  logical  \ariable  'output"  was  added  to  the  input  file  in  two 
modes  "  true"  or  'false"  When  this  variable  is  set  to  'false',  the  screen  output 
is    not  printed  and  the  required  outputs  are  redirected  to  output  tiles 

The  following  list  describes  the  output  files  ana  the  uata  they  contain 
fort. 2      This  is  for  user  supplied  airfoil  coordinates,  if  desired 
fort. 3       This  tile  contains  the  global  coordinates  of  AF1    at  each  time  step. 
fort. 4      This  file  contains  the  global  coordinates  of  AF2  at  each  time  step, 
fort. 7      This  file  contains  the  lift,  moment,  and  drag  coefficients  for  Doth  airfoils 
fort. 8  :    This  file  contains  the  pressure  coeff.  for  AF1   at  each  time  step. 
fort. 9     This  file  contains  the  pressure  coeff  for   AP2  at  each  time  step  . 
fort.  10    This  file  contains  the  first  airfoil's  core  vortex  (wake)  positions, 
fort.  11    This  file  contains   the   second    airfoil's  core  vortex  (wake)  positions, 
fort.  12:  This  file  contains  the  angle  of  attack  at  each  time  step, 
fort.  13.    This  file  contains  the  y-axis  translational  motion  at  each  time  step. 

C.    LOEWY'S  ANALYSIS 

Loewy  [4]  developed  a  two  dimensional  model  for  the  unsteady  aerodynamics 
of  the  blades  of  a  hovering  rotor,  including  the  effect  of  the  shed  wakes.  Figure 
(4  1).    from  reference  [40],    shows    his    two-dimensional  model    of   the    helical 
wakes.    The  wake  in  the  case  of   a    multi-bladed  rotor  in    hover     consists  of 
helical  vortex  sheets  below  the  disk,  one  from  each  blade. 
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Assuming  sinusoidal  motion,  the  ciownwasn  over  the  reference  airfoil.  V,  and 
its  vorticitv  (either  bound  or  shed)  mav  be  written,   respectively  as 


V     =    v\e;"1     .    v    =  v  e   u  i4  1) 


where  the  barrea  quantities  are  complex  The  induced  velocity  at  a  point  ;  on 
the  airfoil  resulting  from  an  element  ot"  vorticitv  of  strength.  ;10  at  a  general 
point  in  the  wake  is  given  by  [4] 

y     (x-c)dc 

'  n  vi  - 

dw(x.t)   =  ;--  ;  (4.2) 

2*  {  (x-c)~-  (nQ-rq)-(h)-] 

where  O  is  the  total  number  of  blades  and  h  is  the  vertical  distance  between 
successive  rows  of  vorticitv,  n  is  the  number  of  revolution  index,  and  q  is  the 
number  of  blades  index  as  shown   in  Figure  4  1.    The  wake  spacing,  h  is  defined 
as  [4]: 


27IU 


(4.3) 
bQQ 

where  u    is  the    inflow    velocitv    (  time-averaued  normal  to  the  disk). 
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The  totaJ  downwash  over  the  airfoil  can  be  expressed  conveniently  by    summing 
the  integrals  involving  the  bound   voracity   and  the  row  of  vorticity    in   the   plane 
of  the  reference  blade  with  those  for  the  rows  of  vortices  below   the  plane  of  the 
rotor 
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The  first  two  terms  on  the  right-hand  side  of  this  integral  equation  are  those 
that  arise  in  the  classical  unsteady  airfoil  theory  of  Theodorsen  [  3  1.  the  third 
and  fourth  terms  contribute  the  influence  01  all  the  voracity  below  the  Diane 
ot  the  rotor  disc.  Loewv  solved  the  integral  equation  using  Sohngen  s  inversion 
formula  [31].  and  obtained  his  lift  deficiency  function.  C  (k.m.h).  which  is  given 
bv 

C  (k.m.h)  -  F  (k.m.h)  -  i    G  (k.m.h)  (4  5; 

where  k.  is  the  reduced  frequency  (defined  bv  )b/l).  m  is  given  bv  the 
ratio  co/O,  h  is  the  wake  spacing  defined  before,  and  the  real  ana  imaginary 
parts  oi   Loewy's  lift  deficiency  function.  F  .  G    are  given  by  [4] 

J.ctA-  (Y,  -(3J.)B 


F'(k.  m.  h) 


A"  -  B" 


i 


(Y    rpJ.)ArJaB 
-G' (k.m.h)   = 

A"  -  B" 

where 

A  =  J,a  +  Y+JB  ,      -B  =  Y.  -i-J.p-J  a 

1  o        ir  1  r  o 

and 

kh        -kh  _    .     /_ 

e      -  e  „  2sin  (_~m ) 


a  = .3 

kh      ^  .  _.         .  -kh  kh      -,  . -,         ,  -kh 

e     -  2cos  (2:im)  +  e  e      -  _cos  (2n:m)  +  e 
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Here  it  can  lie  seen  that  as  the  spacing  between  rows  of  vorticitv,   h, 
becomes    infinite.   F'   approaches  F  and  G'   approaches  Ci.   Also  as  k    tends  to 
infinity,    C'(k.m.h)     approaches  C(k).   When  h  is  zero,  all   the  vorticitv  lies  in 
the  plane  of  the  airfoil,    singularities  result,  and  the  meaning  of  anv    obtained 
results,  phvsicallv  speaking,  is  not  clear  [4]      Loewv  [4],  and  Couch  [40]  showed 
plots  for  F    and  (i      versus    reduced    frequency  k    Loewv    treats  an  infinite 
number  of  wakes  and  Couch  explores  a  finite  number  including  the  special  case 
of  a  single  wake.     These  plots    show  the  effects    due  to  wake  phasing,  m  ,  and 
wake  spacinu,  h.    A  case  of    special  interest  is  the  wake  sheddirm  from  one 
preceding  blade  only,    lor  this  case  Loewy's  lift  deficicnev  function.  C  .  reduces 
to  the  following  lift     function.  C*,    obtained    by     Wood  and  Couch  [40] 


C*    =   F*    -  i  G  *  (4  6) 


where        F     =  — : 


^,V+#V 


G,  =    ,(V^W^-/|(I-2^V)^ 


/7=l 


The  propulsive  force  coefficient  can  be  calculated  and  the  results,  obtained  by 
Wood  and    Couch  [40],  using  F*  and  G*  are  : 


For  pure  plunge   .      Cp    =  77 Jr n'^F')2  +  (C7'): 

For  pure  p,tch  ^  _^:  ^ )2^.)2)m_  ^  _#)]     ^  _ 


[k2 


2  k11    [2       1    k 
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D.  EVALUATION  OF  THE  TWO-FOIL  PANEL  CODE 

Several  analyses  are  available  in  the  aeroavnamic  i'terature  which  permit  a 
: miteci  evaluation  of  the  two- foil  panel  code.    These  are: 

^round    Effect  on  Airfoil  Lift  in  Steady   Flight 

b.  ...i,  ika  et  al  [44]  analyzed  the  ground  effect  on  a  2-D  flat  piate  using 
the  method  of  conformal  transformations  Figure  -1  2  a  shows  the  panel  code 
results.  Figure  4.2  b  reproduces  Tomotika  results,  h  seen  that  the  two  methods 
are  in  good  agreement.  Figure4  3a  and  Figure  43b  show  the  comparison  be- 
tween the  result  obtained  by  vVegley  [46]  for  a  i0%  thick  RAE101  airfoil  in 
ground  effect  and  the  panel  code  prediction  for  a  NACA  0010  at  zero  angle  of 
attack.  The  two  results  are  again  seen  to  be  in  good  agreement  considering  the 
fact  that  the  RAE101  airfoil  (whose  coordinates  are  not  known  )  was  approxi- 
mated by  a  NACA  0010  airfoil. 

2.  Bip.^.e  Effect 

Glauert  [45]  defined  a  factor  B  which  gives  the  lift  reduction  on  the  airfoil 
in  an  unstaggered  biplane  configuration  compared  to  the  lift  on  a  monoplane  at 
the  same  incidence.  Table  4.1  gives  the  correction  factor  for   several  non-dimen- 
sional biplane  spacings.  It  is  seen   from  this  table  and  from  Figure  (4  4)  that  the 
two  results  are  anain  in  reasonable  aureement. 
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Figure  4.2  a  CI/CIo  vs  h/c  (   Panel   (Ode) 
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Figure  (4.2  I))  CI/CIo  vs   h/c  (    From    Rcf.   [44|  ) 
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Figure  (4.3a)    Local  Velocity   Distribution,  V(/V  vs  X/C     (Panel  Code) 
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Figure  (  4.  3  b)    Local   Velocity   Distribution,  V(/V  vs  X/C    (after  Wegley  |46|) 
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Figure  (4.4)    Collection   Factor,   B  vs  the  Chord-gap  Ratio  h/c 
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TABLE  4.1     CORRECTION   FACTOR   FOR  A   BIPLANE 


li/c 


B  (  Code  )         B  (  rcf.   [451  ) 


1 

0.50 

r 

0.770 

0.730 

! 

1 
i 

0.57 

0.847 

0.800 

1.00 

0.902 

0.855 

1.25 

0.942 

0.895 

! 

1.50 

0.97 

0.920 

.       _ 

I 

Another  ground   interference  study  is  due  to  Bagley  [46].  As  can  be  seen  from 
the  comparison  of  the  panel  code  (Figure  4.5  a)  with  his  predictions  (  Figure  4.5 
b  ),  in  ground  effect,    the  lift  is  reduced  over  the  upper  surface  and  increased 
over  the  lower  surface     As  the  airfoil  approaches  the  ground,  the  (low  between 
the  airfoil  and  ground  is  reduced  and  hence  the  pressure  is  increased. 
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Cp  vs  X/C  (SS.h/c- 5.  A0A=1  1  cleg.) 


i 


Figure  (4.5  a)    Pressure  Distribution  vs  X/C    (Panel  Code) 


1  II           1 

ll            .  ... 

r,inr.i  f   mi   ii  Air 

p 
1  'j 

ll 

1  1 
1  \ 
\\ 

\   \ 

M  A  1     I'l  A  IE    tlf  AH 
ciif  mill) 

Mllllll    All      ACIIOf  Oil 
<i 

1  ( < 

\       \ 

»-/r  -1/2                      Vh 
'<  =||* 

M   '• 

unit  it 

S^^       -             C 

sum  ACE 

-~^^  „ 

r  , 

'           '            ^"""^S 

O  7 
HIKE  It 

ni         o  <s       0_8.  ,  - -J 

O   S 

surface. 

1   fl 

_j i      .  .     i 

Fi8'"-c  (4-5  b)    Pressure  Distribution  vs  X/C    (   Baglev 


y) 


73 


E.  ANALYSIS  OF  OSCILLATORY   BLADE   INTERFERENCE   EFFECTS 

Having  determined  in  Chapter  III  that  pure  flapping  of  a  single  helicopter 
rotor    blade  is  unlikely    to  produce  a  drag  reduction    sufficient  to    explain    the 
OH-6A    HHC  power  reductions  [  1  ].  solely,  due  to  the  Katzmayr  effect,  we 
investigate  in    this  section  the  possibility    ot'   favorable    blade/wake  interference 
effects.    Consider  the  arrangement  of  the  two  airfoils  shown  m  Figure  (4  6)    The 
wake  shed  from  airfoil  -2    impinges   airfoil  w  1.    if  the  vertical  spacing  between 
the  two  airfoils    falls    below  a  certain  value,  depending  on  the  amniituue  of 
oscillation  ot  airfoil  #  2  .  This  case  cannot  be  analyzed  with  the  unsteaav  panel 
code  because  the  vortices  start  to  penetrate  the  second  airfoil  or  because  the 
wakes    shed    from  the  two  airfoils  come  in  contact  with  each  other,  as  shown 
in  Figure  (4.7). 

On  the  other  hand,  if  the  vertical  non-dimensional  spacing  is  two    or  more 
no  problems  occur.    A  second  important  parameter  is  the  horizontal    spacing  be- 
tween the  two  airfoils.    This    parameter   controls    the  phasing  between    the  two 
wakes.    Hence  the  unsteady  two-foil  code  permits  the  analysis  ot    Loewys  ar- 
rangement for  the  special  case  of  two  wakes.    This  case  is  considered  by  Wood 
and  Couch  [40]  in  closed  form  and  presented  in  Ref   [39].    Furthermore,  the 
single  airfoil  case  is    recovered    if   the  vertical  spacing  is  chosen  to  be    large 
because  the  effect  of  the  first    airfoil  on   the  second  must  diminish  with  increas- 
ing   distance  between  the  two  foils. 
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Figure  (4.6)     Wake  Pattern  (  Airfoil  #  2   Plunging  at   High  reduced   Frequency 
and  Close  (o    Airfoil  #  I  ) 
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Figure  (4.7)    Wake   Interference  Between  Two  Plunging  Airfoils  in  Close 
Proximity 
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Figure  (4  8).  and  Figure  (4.9)  show  the  computed  lift  and  drag  (thrust)  co- 
efficient for  airfoil  #  1  in  the  presence  of  airfoil  -  2.  Both  airfoils  are  oscillating 
in  plunge  with  an  amplitude  of  0  14  and  a  reduced  frequency  ot  0.0617  The 
vertical  spacing  between  the  two  airfoils  is  200  In  this  case  the  lift  and  thrust 
coefficients  are  found  to  be  identical  to  the  ones  obtained  from  the  single  airfoil 
code,  thus  showing  that  a  distance  of  200  is  sufficient  to  recover  single  airfoil 
results. 

In  Figures  (4  10)  to  (4  21).  on  the  other  hand,  the  vertical  spacing  is  reduced 
to  2.  For  these  calculations,  the  time  step  is  3  to  .5  time  units  The  plunging 
amplitude  is  again  0  14.  the  reduced  frequency  is  0  0617  The  twelve  figures 
4  10  to  4  21  show  four  different  cases  of  phasing  between  the  two  wakes.  In 
Figures  (4.10)  to  (4  12),  the  wake  pattern,  lift  and  drag  are  shown  for  zero 
phasing,  m  =  0,  Figures  (4.  13  )  to  (4  15)  show  a  phase  of  m  =  0.20833,  Figures 
(4.16)  to  (4.18)  show  the  case  of  m  =  0.25,  and  Figures  (4  19)  to  (4  21)  show 
the  case  of  m  =  0.5  . 

A  closer  inspection  of  these  figures  reveals  the  strong  effect  ot  the  wake 
phasing.  For  m  =  0,  airfoil  #  1  experiences  a  reduction  in  thrust  compared  to  the 
single  foil,  whereas  for  m  =  0  5  airfoil  #  1  experiences  a  significant  increase  in 
thrust.  These  results  are  summarized  in  Figure  (4.22)  where  the  single  airfoil 
results  are  also  shown  as  well  as  the  predictions  bv  Wood  and  Couch  [40] 
when  specialized  to  the  case  of  two  interfering  wakes  only.  The  agreement  be- 
tween the  two  predictions  is  encouraging. 
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Figure  (4.8)    (I  vs  T*   (  plunge,  ...  =  0,  h  =  200  ) 
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Figure  (4.9)    C(,  vs  T*  (  pIunge?  „,  =  0,  h  =  200  ) 
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Pigure  (4.10)     \V;ike   Pattern  (   plunge,  in  =  0,   h  =  2) 
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Figure  (4.11)     CI  vs  T*    (  plunge,  m  =  0,  h  =  2  ) 
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F»8»to  (4A2)       Cd  vs  ,p    (        m        m  =  0    fc 
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Figure  (4.13)    Wake  pattern  (  plunge,  m  =  0.20833,  h  =  2  ) 
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Figure  (4.14)      CI  vs  T*  (plunge,  m  =  0.20833,  h  =  2) 
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Figure  (4.15)      Cd  vs  T*  (  m=  0.20833  ,  h  =  2) 
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Figure  (4.16)   Wake  pattern  (  plunge,  m  =  0.25  ,  h  -  2  ) 
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Figure  (4.17)       CI  vs  T*    (plunge,  m  =  0.25  ,  h  =2  ) 
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Figure  (4.18)     Cd  vs  T*    (  plunge,  m  =  0.25  ,  h  =2  ) 
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Figure  (4.19)      Wake  pattern  (plunge,  m  =  0.5  ,  h  =2  ) 
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Figure  (4.20)     CI  vs  T*  (  plunge,  m  =  0.5  ,  li  =2  ) 
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Figure  (4.21)      Cd  vs  T*    (  plunge,  m  =  0.5  ,  h  =  2  ) 
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Figure  (4.22)    Average  Propulsive  Foree  Due  to  Plunge  Versus  Wake  Phasing,  in 


92 


Similar  results  are  obtained  for  the  pitch  case.    The  airfoils  studied  are  again 
NACA  0007  foils,  pitching  about  the  quarter-chord  point  with  an  amplitude  of 
I   degree  and  a  reduced  frequency  of  0  0617.    Figures  (4.23)  to  (4  32)  display 
the  same  phasing  values  as  in  the  plunge  case.    These  figures  and  the  summa- 
rizing Figure  (4.33)  show  again  that  values  of  m  near    zero    and    one    produce 
increased    drag  compared  to  the  single  airfoil  case,  but  values  of  m  near  0.5 
produce  thrust,  whereas  the  single  airfoil  produces  drag.  Both  analyses,  the  anal- 
ysis by  Wood  and  Couch  [40]  for  the  special  case  of  two  wakes  only    and  the 
panel  code,  produce  similar  trends.   However,  the  quantitative  agreement  is    not 
as  good  as  in  the  plunge  case    The  precise  reason  is  not  sufficiently   understood 
at  this  time     It  appears  that  the  resolution    of  the  suction  peaks  over  the  leading 
edge  and  therefore  the  computation  of  the  thrust  is  more  sensitive  in  the  pitch 
case  than  in  the  plunge  case.    Furthermore,    it  must  be  remembered  that  Loewy's 
analysis  for  an  infinite  number  of  wakes  and  that  by    Wood  and  Couch  [40]  for 
two  wakes  are  based  on  the  flat-plate  assumption,  whereas  the  panel  code  requires 
a  minimum  blade  thickness  in  order  to  obtain  accurate  results.  A  further  evalu- 
ation of  the  panel  code  can  be  done  for  the  drag  produced  by  a  single  airfoil 
pitching  about  the  leading  edge.    In  this  case  the  analytical  theory  of  Garrick  [5] 
and  the  computational  results  of  Bosch  [  15  ]  are  available.    It  can  be  seen  from 
Figure  (4.34)  that  the  trends  are  again  in  agreement,  but  the  code  results  deviate 
from  (jarrick's  and  Bosch's  results  with  increasing  frequency     These  deviations 
must  be  attributed  to  the  geometry  differences  (  flat  plate  versus  NACA  0007  ) 
and  the  higher  amplitude  used  in  the  panel  code.   Note  that  the  pitching  airfoil 
starts  to  develop  thrust  only  at  a  reduced  frequency  greater  than  about   1.8. 
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Figure  (  4.23  )  Wake  pattern  (    pitch,  m  =  0  ,  h  =2  ) 
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Figure  (4.24)       CI  vs  T*  (    pitch,  m  =  0  ,  h  =2 ) 
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figure  (4.25)      Cd  vs  T*  (    pitch,  m  =  0  ,  h  =2  ) 
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Figure  (4.26)     CI  vs  T*  (    pitch,  m  =  0.20833  ,  li  =2  ) 
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Figure  (4.27)    Cd  vs  T*  (    pitch,  m  =  0.2083  ,  h  =2  ) 
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Figure  (4.28)    CI  vs  T*    (    pitch,  m  =  0.25  ,  h  =2  ) 
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Figure  (4.29)      Cd  vs  T*    (  pitch,  m  =  0.25  ,  h  =  2  ) 
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Figure  (4.30)  Wake  pallcrn  (    pitch,  m  =  0.5 ,  li -2  ) 
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Figure  (4.31)     CI  vs    I*    (   pilch,  m  =  0.5  ,  h  =2  ) 
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Figure  (4.32)    Cd  vs  T*    (   pilch,  m  =  0.5  ,  h  =2 ) 
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Figure  (4.33)     Average  drag  vs  wake  phasing  (  pilch  only  ) 
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Figure  (4.34)  Average  drag  versus  reduced  frequency  (  single  airfoil  ) 
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V.  PROPULSIVE  EFFICIENCY 
AND  FLUTTER  ANALYSIS 

A.     INTRODUCTION 

In  the  previous  two  chapters  it  was  shown  that  a  single  airfoil  which 
pertorms    harmonic  plunge  oscillations  generates  a  certain  amount  of  forward 
thrust  (  Katzmayr  effect )  and  that  this  thrust  may  be  enhanced  by  the  oscillatory 
wake  shed  from  the  preceding  blade  if  the  frequency  ratio  or   phasing  m  is    in 
the   range  0.2   <   m   >  0.7   approximately   0.5.    It  was  also  shown  that  a  similar 
enhancement  occurs  for  pitching  blades.  However,  pitching  airfoils  develop  thrust 
only    for   relatively    high   values  of  reduced    frequency,    k   >   0.6   by    Garrick 
results  [5]. 

Furthermore,  it  is  well  known  that  a  pitching  airfoil  may  develop  negative 
aerodynamic  damping  and  thus  experience  single-degree-of-freedom  flutter.    As 
nown  by  Loewy  [4],  the  possibility  of  flutter  is  enhanced  by  wake  interference. 
This  is  an  effect  known    as    wake-induced  flutter.    It  is  therefore  of  interest  to 
study  the  effect  of  various  parameters  on    propulsive  efficiency  and  on  flutter 
instability  in  more  detail. 
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B.  PROPULSIVE  EFFICIENCY 


The  propulsive  efficiency    is  defined  as  the    ratio  of  the  average  work    of 
propulsion  to  the  work  required  to  maintain  the  oscillation.  If  this  factor  is  zero 
or  negative,  we  expect  energy  to  be  extracted  from  the  air  and  fed  to  the  oscil- 
lating system.  For  a  flat  plate  performing  a  flapping  motion  given  bv   hit),  where 
h(t)  =  h0  elwt    the  result  for  the  horizontal  force  per  unit  span,    averaged  over 
one  cycle,  is  given  bv    [42]  : 


•)     ->  r      -i 


D  =  -TrpboTlT  _F~(k)  +G"(k) 


(5.1) 


where  F(k)  and  G(k)  are  the  real  and  imaginary  parts  of  the  Theodorsen  lift 
deficiency  function  C(k).  For  the  case  of  pitching  motion,  the  average  horizontal 
force  per  unit  span  length  is  given  by  the  following  equation,  (5.2): 


D  =  -rcpb  (a~a"\  [F"  +  G  j 


hi- 


"1 

+ 

-F 

_2 

_ 

1  >       F_GH 

2  ~   >      H2  ~  K  ^2 


+  a 


where  ac  is  pitch  amplitude  and  (a)  is  the  position  ot^  axis  of  rotation.  a=  -  1 
corresponds  to  the  leading  edge,  a  =  +  1  is  the  trailing  edge.  The  propulsive 
efficiency  then  is  given  by  : 

DU 


n  = 


w 


(5.3) 


W    is  the  average  work  required  to  maintain  the  oscillation.    For  pitch,  it  is 
given  by    [42]: 


3     ■>        2     1  (\ 

W  =  7ipb  co  UcT  {-  I  —  a 
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a  +  -; 


o       J      k  I 


(5.4) 
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loi   plunge,  it  is  given  by 

W  -  Tupb2!!2^-) %    F(k)    -  Tipbh\)2UF(k)  (5.5) 

Then  the  efficiency  factor,   t],  can    be  expressed  as.  for  plunge  : 

l;   (k)   i  (T(k) 
"'.=  F(k,  (5"6) 

■<uk\    for    pitch 


in  ;      I 


I  '  -a)-   f'   -Gfl  +a 

2       \i        /  |  2  v  2         ■  2       k  ^  2 


(5.7) 


Figure  (5.1)  shows  the  propulsive  efficiency  of  a  single  plunging  airfoil  as 
a  function  of  reduced  frequency  Slowly  flapping  airfoils  reach  very  high  values 
of  efficiency,  approaching  one.  As  the  frequency  of  flapping  is  increased,  the 
efficiency  falls  to  50  percent  Pitching  airfoils,  on  the  other  hand,  start  out  with 
an  efficiency  of  no  more  than  50  percent  at  high  frequency  of  oscillation  which 
quickly  falls  to  zero  and  to  negative  values  with  decreasing  frequency.  A  negative 
propulsive  efficiency  means  that  no  work  is  required  to  maintain  the  oscillation, 
but  instead  energy  is  extracted  from  the  airst ream.  This  implies  a  loss  of  damping 
and  thus  the  possibility  of  single-degree-of-freedom   flutter  in  pitch. 
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Figure  (5.1)     Propulsive    Efficiency     n   versus    1/k   for  a    plunging    airfoil 


109 


Figure  (5.2)  shows  the  propulsive  efficiency  of  an  airfoil  which  pitches  about 
the  leading  edge.    The  efficiency    quickly    drops  to  negative  values    as    the  fre- 
quency is  reduced.    As  seen  in  Figure  (5.2),  infinitely    large  negative   values  are 
reached  as  the  inverse  frequency  approaches  24,  corresponding  to  a  reduced  fre- 
quency   of  0038.    This  is  the  flutter  condition,  as  will  be  explained  in  more 
detail  in  the  next  section. 

Similar  trends  are  found  for  other  pitch  axis  locations,  as  seen  in  Figure 
5  3  (  pitch  axis  is  at  the  quarter  chord  and  mid-chord  point  ).  The  three  quarter 
chord  point  is  an  exception.  Here  the  propulsive  efficiency  is  minus  one  at  very 
high  reduced  frequency,  which  falls  to  larger  negative  values  with  decreasing 
frequency  (  Figure  5.4  ).  However,  it  is  interesting  to  note  that  a  small  shift  of 
the  pitch  axis  of  only  0  5  %  in  either  direction  restores  the  efficiency  to  0.5 
for  large  frequencies,  from  which  it  starts  to  drop  with  decreasing  frequency. 

These  results  show  that  there  is  a  fundamental  difference  between  airfoil 
plunge  and  pitch  motions.  Plunge  oscillations   of  a  single  airfoil  always  produce 
thrust,  and  the  efficiency    decreases  as  the  frequency    is    increased.    Pitch    oscil- 
lations,   on  the  other  hand,  produce  thrust  only    at  relatively  high  values  of  re- 
duced frequency,  k  >  0.6.  This  possibility  of  si ngle-degree-of- freedom  in  pitch 
flutter  is  discussed  in  more  detail  in  the  next  section. 
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ligiiro  (5.2)    Propulsive  Efficiency    i]  versus   l/k  for  an  airfoil      pitching 
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Figure  {5.3)    Propulsive  Efficiency  r\  versus    1/k  for  an  airfoil    pitching 
about  the  quarter  -  chord  point  and  mid -chord  point 
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Figure  (5.4)    Propulsive  Efficiency  n,  versus   1/k    for  an  airfoil 
pitching    about  the  three  -  quarter- chord  point 
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C.  SINGLE-DEGREE-OF-FREEDOM  AIRFOIL  FLUTTER 

Consider  an  airfoil  hinged  at  its  leading  edge  but  elastically  restrained  from 
rotating  about  this  axis  by  a  torsion  spring  with  constant  Ka  (ft-lb/rad).  The 
airfoil  is  placed  in  a  low-speed    airstream  so  that  the  unstrained  position  of  the 
spring  corresponds  to  zero  angle  of  attack    a.    Figure  (5.5).  The  equation  of 
motion  for  this  single-degree-of-freedom  system  is  given  by  : 


I  a  +  K  a  =  M  .  (5.8) 

ix  a  v 


Mv  is  the  aerodynamic  moment  due  to  a  (t)  and  Ia  is  the  moment  of  inertia 
about  the  leading  edge.  Assuming  that  a(t)=  aoe1C0t,  then  equation  (5.8)  can 
be  written  as 

I 

5L    i- 

V  co 


7ipb 


+  m     =0.  (S-?) 
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where  coa    is  the  natural  frequency  of  torsional  vibration  and  is  given  by 
K 


a 


co  (5.10) 

06  M    I 

ll     a 

and  mv  is  shorthand  for  the  dimensionless  aerodynamic  coefficient  and  can  be 
written  as  [41]: 

m    =-        — •- =  M    -(L    +-]fi+al+L.  f-  +  al  .  (5.11) 

>'  .42         icot  a       \    a      -jj  \i        J         h  V?        / 

Tipb  co  a  e  *" 

r  o 
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where "  a  "  is    !he  clastic  axis  position    and    at    the    leading  edge  a  = 


m, 


is  a  complex  number  and  it  is  a  function  only  of  reduced  frequency,   k  =  wb/U. 


Equation  (5.9)  can    he    split  into  real  and  imaginary  parts. 


Re {m    }  -  Re 


.  irpb 


It 


V  (i)   J 


(5.12  a) 


Im   (mv}=0     (5.12b) 

Flutter   occurs  at  that  value  of  the  reduced    frequency    where  the  imaginary 
(  out-of-phase)   part  of  the  aerodynamic  moment  becomes  zero,  provided  that  the 
corresponding  real  (in-phase)  part  yields  a  non-imaginary  flutter  frequency. 


Figure  (5.5)  An  airfoil   restrained  to  rotate  about  its  leading  edge 


in  two-dimensional   flow. 
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Figure  5.6  shows  the  variation  of  the  real  and  imaginary  (  damping  )  parts  of 
the  pitching  moment  for  an  airfoil  pitching  about  the  leading  edge.  It  is  seen 
that  the  pitch  damping  goes  through  zero  at  a  reduced  frequency  of  0  038  This 
value  agrees  with  the  previously  found  value  using  propulsive  efficiency  ( energy 
extraction  )  considerations. 

Figure  5.7    displays  the  comparison  between  the  pitch  damping  coefficients 

(note  that  the  pitch  damping  coefficient  is    defined    by    Theodorsen  [3]    as  : 
2  ~                                                                      2 
C     =  ~—    M    -fi+a](L    +M.)  +  f-+a]  L.     (5.13) 

m  2  a       v?         /        ex  lr        v.?         /hi 

The  difference  in  Cm  shown  in  Figure  4.8  is  due  to  a  difference    between   the 
Ma  ,  La ,  Mh   and   Lh   terms   in    Equation  5.13    as   determined   by    panel  code 
(  Riester  [43]  )  and  Theodorsen's  theory  [3].  The  imaginary   part   in  Figure  5.7 
is  obtained    by    Equation  5.11  computed  by  the  panel  code  and  Theodorsen's 
analysis.    It  is  seen    that  the  two    computations    are  in  reasonable  agreement, 
recognizing  the    difference    in    geometry    (  NACA  0007  versus  flat  plate  ). 

Systematic  variation  of  the  pitch  axis  location  shows  that  zero  pitch  damping 
is  possible  for  axis  locations  upstream  of  the  quarter -chord  point.  No  flutter  is 
possible,  on  the  other  hand,  for  axis  locations  at  or  downstream  of  the  quarter 
chord  point,  Figure  5.8.  Loewy  [4]  has  shown  that  wake  interference  greatly 
increases  the  possibility  of  flutter.  The  unsteady  panel  code  provides  an  oppor- 
tunity to  investigate  and  to  extend  this  well-known  finding.  This  is  a  problem 
which  should  be  studied  in  more  detail. 
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Real  6   Imag.  pirt  of  m     vs  k  at   x  =  -  1 
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I'iRiirc  (5.6)  Variation  with  reduced  frequency  k  of  (lie  real  and  imaginary 
parts  of  the  non-dimensional  aerodynamic  moment  ui^  due  to  pitching  of 
an  nil  foil  about  its  leading  edge  in  an  incompressible  flow. 
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Cm  vs  k    ( pitching  airfoil,   aQ   =   1°) 
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Fipurc  (5.7)    Comparison  of  Cm  (Imag.)  vs    Kp  for  NACA0007  Oscillating 
in   pitch  about  (he  leading  edge  for  a0  =  1°   and  IN  =200  panels. 
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Figure  (5.8)  Imag-  (my )  vs  k  (for  an  airfoil  pitching  at   ca  =  -0.5,0.0,0.5,  and  1.0) 
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VI.    RESULTS  AND  CONCLUSIONS 

A  DISCUSSION  OF  RESULTS 

This  research  was  undertaken,  to  study    the  effect  of  unsteady  aerodynamics 
on    reductions    of   main    rotor  shaft  torque  and    engine  power    measured    when 
the  HHC  system  was  applied  to  the  OH-6A  rotor.   In  Chapter  III  we  study  pro- 
pulsive force  resulting  from  pure  plunge  oscillation  of   an  airfoil  in  in  viscid, 
incompressible  flow  The  results  of  flat  plate  theory,  Theodorsen  [3]  and  Garrick 
[5],  are  used  to  validate   panel  code  results.    Comparison   with  Garrick   showed 
good  agreement  in  pure  plunge,  especially,  for  lower  values  of  reduced  frequency 
and  amplitudes  of  oscillation. 

Figures  3.5  to  3  ,    show    the  average  drag  coefficient  versus    plunging 
amplitude  for  three  values  of  reduced  frequency,  k=  0.1,  1.  and  2  respectively. 
It   is    seen    from    Figure  3.5,    (  reduced    frequency,  k    is   0.1  )   that  agreement 
between    the  panel  code  and  Garrick  [5]  is  very    good    for  a  range  of   non- 
dimensional  plunging  amplitude  up  to  15  to  20  %  of  the  blade  chord.     On  the 
other  hand,  for   higher   values  of  k,  the  two  results  are  in  good  agreement  for 
a  range  of  plunging  amplitude  up  to  approximately  10%.  Single  airfoil  code,  was 
used  to  compute  the  time  history  of  the  drag  (thrust)  coefficient  of  the  NACA 
0015  airfoil  (the  blade  section  of  the  OH-6A  ). 
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Figure  3.12  and  Table 3. 3   show  the  average  drag  coefficient  versus    reduced 
frequencv    for  an  NACA  0015  airfoil  which  executes  plunging  oscillation  with 

1  "  plunging   amplitude.    The  OH-6A  data    presented    in  Table  3.1,  Wood    [42], 
and    Figures  3.8    to    Figure  3.11  from    Ref.  [27],  were  used    to    estimate  the 
propulsive  force  per  blade  obtained  due  to  plunging  oscillation  for  1  P  frequency 
case.     In  these  figures,  the  difference   between    the  deflections    when  the  HHC 
is  "  on  "  and  "  off  "  represents  the  plunging  amplitude.  The  controller  phase  is 
also  significant.  Phase  angle  of    90°  or    120°  are  those  where  helicopter  vibra- 
tions were  most  sever,  whereas  at  phase  angle  of  300°  or  330°.  the  helicopter 
vibrations  were  significantly  reduced. 

We  look  for  the  greatest  power  benefit  to  occur  at  the  phase  angle  of 
lowest  vibration,  that  is  300°  or  330°.  Table  3.2  shows  the  amount  of  propulsive 
force  and  torque  generated  at  each  blade  segment.  The  results  of  this  table 
showed  that  the  contribution  of  IP  (first  harmonic)  were  3  horsepower.  Wood 
et  al  [27],  applied  the  Garrick  equation  and  considered  the  effect  of  blade  plunge 
response  in  relation  to  the  second  to  twelfth  harmonic  loads.  He  found  for  2  P 
harmonic  only  (HHC  "on",  330°  phase  at  60  knots  ),  the  power  gained  was 
about  11.3  horse  power.  For  more  detail,  the  reader  is  encouraged  to  study  Table 

2  of   Ref.  [27]  l  . 


Wood  E.  R.  .  Iligman.  J.  and  Ramesh  Kolar.  "  Higher  Harmonic  Control 
Promises  Improved  Dynamic  Interlace  Operations"  AGARD  Proceedings  of 
the  78  th  Flight  Mechanics  Panel.  May    1991.  Seville.  Spain. 
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Chapter  IV,  considers  the  influence  of  the  rotor  wake  interaction  and  examines 
whereas  such  interaction  might  cause  favorable  effects.  First,  the  case  of  an   iso- 
lated   airfoil  (  single  airfoil  oscillating  in  plunge)  with  absence  of   wake  was 
considered  ( or  the  wake  located  at  great  distance  from  rotor  ).  This  is  the  "Katz- 
mayr  effect"  case  where  the  resulting  drag  force  is  a  propulsive  force.    It  was 
the  primary    case  considered    by    Garrick  [5]. 

The  advantage  of  the  Panel  Code  is  that  it  provides  us  with  time  histories 
of  unsteady  lift,  drag,  and  the  wake  trajectory  bv  which  we  can  better  understand 
and  analyze  the  problem.  In  each  case,  the  wake  pattern  was  plotted  in  order 
to  be  sure  that  the  wake  produced  by  the  second  airfoil  would  not  interfere  with 
the  reference  airfoil.  Also  the  lift  and  drag  (or  the  propulsive  force)  were  plotted 
versus    non-dimensional  time,  T*  (  T*  =  t.U/c  ). 

The  closed  form  solutions  are  based  upon  a  very  thin  flat  plate  airfoil 
undergoing  very  small  motions  when  we  model  it  in  the  unsteady   panel   code. 
For  the  calculations  an  NACA0007  airfoil  was  selected.    The  airfoil  sets  at  zero 
mean  angle  of  attack  and  the  reduced  frequency,  k,    is    =  0.0617 

For  the  OH-6A,  the  value  of  reduced  frequency  ,  k  at  70  %  radial  station , 
was  calculated  and  found  to  be  .1234  for  the  32  Hz  exciting  frequency.  The 
reference  length  used  by  the  panel  code  is  the  blade  chord  or  2  b,  so  the  0.1234 
reduced  frequency  corresponds  to  0.0617  per  the  panel  code  definition.  Non- 
dimensional  wake  spacing,  for  single  airfoil  analysis,  was  h  =  200  (  mathemat- 
ically equivalent  to  infinite  spacing,  or  we  say  that  the  interaction  between  the 
two  airfoil  is  very  small)  . 
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The  amplitude  of  plunge  oscillation  was  h()  =  0.14     (  by    the  panel  code  = 
.07,  tor  the  same  reason  ).  The  panel  code  results  are  compared  to  other  refer- 
ences  in  order  to  know  the  code  limits.  The  comparison  of  the  code  results  to 
that  obtained    by     Ref.   [44],  is  shown  in  Figure  (4.2).  This  figure  shows  the 
relative  lift  coefficient  of  an  airfoil  (  a  wing  of  relatively    long  span  )  as  it  ap- 
proaches the  ground  for  different  values  of  angles  of  attack.    The  figure,  showed 
that  the  results    estimated  by  the    panel  code  and  the  results  obtained  by    Ref 
[44],    are  very  close.      For    h/c  <  0.5    the  results    obtained  bv  the  panel  code 
are  suspected. 

Another  comparison  is  shown  in  Figure  (4.4).  This  figure    and  Table  4  1 
show  the  correction  factor,  B,  versus  the  chord-gap  ratio  in  a  biplane  configura- 
tion. The  correction  factor,  as  defined  by  Ref.  [45],  is  the  reduction   of  the  lift 
coefficient  of  an  unstaggered  biplane  compared  with    that  of  monoplane  at  the 
same  angle  of    incidence.  On  Figure  (4.3)  and  Figure  (4.5),  the  local  speed  ratio 
and  pressure  coefficient  vs  X/C     are  plotted  and  showed  very  good  agreement 
with    Wegley    [46]  and    Bagley    (1959),  respectively. 

The  panel  code  is  suspect  whenever   the  wake  generated  by  an  airfoil 
impinges  on  another  airfoil.    The  numerical  panel  code  method  also  encounters 
difficulty  with  discontinuities.  For  these  reasons  it  is  important  to  carefully  check 
the  time  history  plots  after  each  computer  run. 
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In  Figure  (4.8),  the  time  history    of   the  lift  coefficient.  C]    variation  is 
plotted.      Shown    in    Figure  (4.9    is    the  drag    coefficient.  Cj     versus  the    non- 
dimensional  time  T*     We  observe  from  Figure  (4.8)  that  the  lift  varies  as  the 
frequency  of  oscillation,  with    a  mean  lift  value  of  zero  as  we  expect.  Amplitude 
of  oscillatory    lift  was  found  to  be±  0.0510.    Figure  (4.9)  clearly  shows  the  " 
Katzmayr  effect".    We  noted  that  the  calculated  drag  varies  at  twice  the  plunge 
oscillation  frequency  of  the  airfoil.    Also,  observe  that  the  airfoil  is  generating  a 
propulsive  force  as  indicated  by  the  mean  value  of  drag  coefficient,  which  is 
seen  to  have  a  negative  value. 

Looking  to  the  computer  printout  at,  at  the  non-dimensional  time  T*  = 
0.0,     steady  state,  the  drag    generated    by    the  first    airfoil    was    found    to    be 
Cdss    =  00013 1  for     NACA0007,  which  is  the  airfoil  type  that  was  used  for 
most  of  our   analyses    at   zero  angle  of  attack ,  and     a  total  number  of  panels 
N=  100  (  i.  e  50  panels  at  the  upper  surface  and  50  panels  at  the  lower  surface). 
For   NACA0009,  (this  airfoil  was  used  to  compare  the  code  results  with  that  of 
Bosch's  for  a  pitching  case),  the  steady  state  drag  Cdss  was  found  to  be  .000196 
at  N=  100  and  .000113  at  number  of   panels  N=150. 

It  is  known  that,  for  an  airfoil  at  zero  angle  of  attack,  in  an  inviscid 
incompressible  flow  field,  steady  state  drag  is  zero,  so  in  our  calculations  that 
steady  state  drag  was  subtracted  from  the  drag  generated  at  each  time  step.  Shown 
in  Figure  (6.1)  is  the  steady  state  drag  coefficient  vs  angle  of  attack  [37].  Figure 
(6.1a)  shows  the  variation  of  steady  state  drag  coefficient  vs  angle  of  attack  for 
three  different  airfoils;  NACA0003,  NACA0007,  and  NACA0012  at  total  number 
of  panel  of  200. 
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In   this  graph,  it  is  shown  that,  for   the  thin  airfoil  NACA0003    at  zero  angle 
of  attack,  the  drag  coefficient  is  very  small  but   as  the  angle  of  attack  increases, 
the  drag  coefficient  highly  increases  more  than  in  the  cases  of  NACA0007  and 
NACA0012.  Figure (6. lb)  shows  the  etfect  of  increasing  the  total  number  of  pan- 
els N  on  the  drag  coefficient  of  the  airfoil     NACA0012.    As  the  total  number 
of  panels  increases,  the  steady  state  drag  coefficient  accumulated  error  decreases. 

Figures  (4.10)  through  Figure  (4  21),  show  the  same  case  as  above,  that 
of  an  airfoil  oscillating  in  plunge,  but  with  the  introduction  of  a  wake  layer 
in  the  near  proximity  of  the  reference  airfoil  (non-dimensionai  wake  spacing,  h 
=2  )  for  wake  phasing,  m  =  0,  0.20833,  0.25,  and  m  =  0.5.  These  correspond 
to  phase  angles  of     0°  ,  75°,    90°  and   180°,  respectively. 

We  see  that  the  phase  angles  relate  to  the  phase  relationship  between  the 
reference  airfoil  and  the  wake  layer  of  shed  vorticity.  These  phase  angles  were 
obtained  in  the  panel  code  by  establishing  the  proper  time  phase  between  the 
initiation  of  the  wake  vorticity  and  the  oscillation  of  the  reference  airfoil.  This 
time  phase  is  in  turn  a  function  of  the  frequency  of  oscillation  which  is  directly 
related  reduced  frequency,  k.  For  the  results  presented  here,  for  k=0617,  this 
corresponds  to  50.92  time  units  for  m=0.0  (0°  phase  );  61.53  time  units  for  m 
=  0.2083  (75°  phase);  and  76.38  time  units  for  m  =  0.50(  180°  phase),  etc. 
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Figure  (6.1)   Drag  coefficient  vs  angle  of  attack 
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If   the  reduced    frequency  .  k,    changes,  this  corresponds  to  a  change  in 
frequency   (at  given  chord  length  and  free  stream  velocity  ).    For  k  =  0  1234.  this 
correspond  50.92  time  units  for  m  =  0.0,  56.23  time  units  for  m  =  0.2083,  and 
63.65  time  units  for  m  =  0.5,  etc. 

Figure  (4.10)  shows  the  wake  position  for  m=  0.0,  the  case  where  the 
wake  shed  from  the  second  airfoil  at  h  =  2  is  in  phase  with  wake  shed  from 
the  reference  airfoil  itself  In  Figure!  4  21  ),  the  time  history  of  lift  variation 
is  shown 

This  is  indicated  by  the  solid  line  for  C\\  .  The  dashed  line,  CU,  indicates 
the  variation  of  lift  on  a  second  airfoil  located  at  the  point  of  initiation  of  the 
wake  layer  at  h  =  2.  Our  discussion  will  focus  airfoil]  ,  indicated  by  the  solid 
line  in  the  figures  mentioned  above.  For  an  airfoil  oscillating  about  zero  angle 
of  attack,  the  mean  lift  is  zero.  We  found  that,  while  the  mean  lift  is  zero,  the 
amplitude  of  lift  oscillation  is  smaller  in  the  presence  of  the  wake  at  zero 
phasing.  From  these  figures,  we  conclude  that  the  effect  of  wake  vorticity  at 
zero  phasing    is  to  reduce  the  "  Katzmayr  effect". 

Comparison  of  oscillatory     lift  values  of  the  75°    wake  phasing    case,  with 
the  case  of   no  wake  or  single  airfoil,    shows  that  they  are  identical,  since  the 
values  at  75  degrees  phase  are  the  same  as  that  for  infinite  wake.  Table  6.1 
proves  this  result. 
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The  time  history  of  drag  is  seen  to  vary  as  twice  the  plunge  oscillation 
frequency  of  the  airfoil.  We  also  observe  that  the  airfoil  is  generating  propulsive 
force  as  indicated  by  the  mean  value  of  drag,  which  is  seen  to  have  a  negative 
value.  Results  of  75  degree  wake  phasing  case,  confirm  that  there  is  a  phase 
angle,  m,  at  which  the  effect  of  the  wake  vorticity  on  the  reference  airfoil  gives 
the  same  results  as  when  the  wake  is  removed  It  has  been  shown  at  a  phase 
angle  of  zero  degrees  that  the  shed  wake  vorticity  increases  the  drag  on  the 
oscillating  refe:ence  airfoil,  i.e.  diminishes  the  "Katzmavr  effect" 

Conversely  we  found  that  at  a  phase  angle  of  about  180  degrees,  there  is  a 
significant  decrease  in  drag  at  the  reference  airfoil  due  to  the  wake  vorticity. 
This  implies  an  intermediate  phase  angle  at  which  the  drag  is  identical  to  the 
"Katzmavr  effect".  It  turns  out  that  such  a  phase  angle  exists.  In  fact,  at  this 
phase  angle  the  lift,  drag  and  pitching  moment  are  found  to  be  identical  to  the 
case  of  no  wake,  and  match  the  "no  wake"  values  identically. 

Table  (6.1)    show  the  aerodynamic  coefficients  Cl(l),  Cl(2),  Cm(l),  Cm(2), 
Cd(l),  and  Cd(2)  vs.  Non-dimensional  time.  This  table  is  a  part  of  one  of  the 
output  files.  It  shows  these  aerodynamic  coefficient  vs   time  from   the  nondimen- 
sional  time  unit  113.3  until   159.1.  The  phase  angle  at  which  this  occurs  is  75 
degree  or  m  =  0.20833,  for  the  case  considered  here. 
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Tabic  (6.1)  Aerodynamic  coefficients  vs  Non-dimensional  time  (from  time 
step  113  to  159)  for  the  case  of  wake  phasing,  m  =75  degree,  wake  spacing 
h  =  2.  reduced  frequency  k  =0.0617  and  plunging  amplitude  hhar  =  0.14) 


TIME       CL( 1 )      CL( 2)      CM(1)      CM(2)      CD(1)     CD ( 2 ) 

113.290932   0.006952   0.0149*5  -0.001387  -0.003432  -0.000011  -0.000039 
114.56386S  -0.001073   0.007130   0.000665  -0.001427   0.000000  -0.000011 
115.836800  -0.009001  -0.000881   0.002701   0.000614  -0.000009   0.000000 
117.109734  -0.016877  -0.008871   0.004674   0.002639  -0.000036  -0.000008 
118.3R2668  -0.024266  -0.0166*4   0.006534   0.004599  -0.000079  -0.000035 
119.655602  -0.031067  -0.024005   0.008236   0.006446  -0.000133  -0.000078 
120.928535  -0.037111  -0.030777   0.009738   0.008134  -0.000194  -0.000132 
122.201469  -0.042247  -0.036791   0.011001   0.009621  -0.000255  -0.000193 
123.474403  -0.046347  -0.041898   0.011994   0.010872  -0.000310  -0.000254 
124.747337  -0.049307  -0.045974   0.012693   0.011855  -0.000355  -0.000309 
126.020271  -0.051051  -0.048918   0.013077   0.012546  -0.000384  -0.000353 
127.293205  -0.051534  -0.050659   0.013140   0.012929  -0.000394  -0.000382 
128.566132  -0.050746  -0.051152   0.012878   0.012994  -0.000386  -0.000393 
129.839066  -0.048701  -0.050385   0.012296   0.012738  -0.000358  -0.000384 
131.112000  -0.045454  -0.048379   0.011412   0.012169  -0.000315  -0.000358 
132.384933  -0.041083  -0.045182   0.010245   0.011301  -0.000260  -0.000315 
133.657867  -0.035698  -0.040872   0.008825   0.010155  -0.000199  -0.000261 
134.930801  -0.029433  -0.035556   0.007188   0.008758  -0.000138  -0.000200 
116.203735  -0.022446  -0.029364   0.005375   0.007146  -0.000082  -0.000139 
137.476669  -0.014910  -0.022450   0.003431   0.005358  -0.000038  -0.000084 
MR. 749603  -0.007014  -0.014982   0.001404   0.003437  -0.000009  -0.000040 
140.022537   0.001046  -0.007147  -0.000655   0.001432   0.000001  -0.000011 
141.295471   0.009072   0.000865  -0.002696  -0.000608  -0.000007   0.000000 
142.568405   0.016865   0.008856  -0.004668  -0.002634  -0.000035  -0.000008 
141.841359   0.024230   0.016627  -0.006521  -0.004593  -0.000078  -0.000035 
145.114273   0.030990   0.023991  -0.008212  -0.006441  -0.000133  -0.000078 
146.387207   0.036981   0.030761  -0.009700  -0.008128  -0.000194  -0.000132 
147.660141   0.042057   0.016775  -0.010948  -0.009616  -0.000255  -0.000193 
148.913075   0.046093   0.041883  -0.011926  -0.010867  -0.000310  -0.000254 
150.206009   0.048996   0.045960  -0.012610  -0.011850  -0.000354  -0.000309 
151.470943   0.050693   0.048903  -0.012985  -0.012541  -0.000382  -0.000353 
152.751877   0.051145   0.050643  -0.013039  -0.012923  -0.000392-0.000382 
154.024811   0.050343   0.051136  -0.012775  -0.012988  -0.000384  -0.000393 
155.297745   0.048306   0.050371  -0.012196  -0.012733  -0.000357  -0.000384 
156.570679   0.045084   0.048364  -0.011319  -0.012164  -0.000314  -0.000357 
157.843613   0.040757   0.045167  -0.010163  -0.011296  -0.000260  -0.000315 
159.116547   0.035427   0.040858  -0.008757  -0.010150  -0.000199  -0.000260 
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The  proper  phasing  is  achieved    by  positioning  the  reference  airfoil  76.38 
time  units  to  the  rig.     of  the  second  airfoil,  the  airfoil  that  generates  the  wake. 
We  found  in  the  case  of  the  plunging  airfoil  that  the  wake  at  zero  degree  phase 
increases  the  drag    force,  while  the  wake  at  75  degree  phase  gives  us  close 
results  to  "no  wake  at  all"  case.  So,  we   expect  that  the  wake  at  180  degree  to 
increase  the  propulsive  force.    Shown  in  Figure  (4.22),  is  the  variation  of  the 
average  propulsive  force  coefficient  versus  wake  phasing,  m.    The  circles  in  the 
graph  represent  panel  code  results  for  the  parameters  mentioned  above. 

Superimposed  on  the  plot  of  the  figure    are  the  analytical  closed  form 
results  of  the  Loewy    theory  as  modified    for  one  wake,  as  shown  in  Figure 
(4.22)   from   Wood  et  al  [38].     Note  the  good  agreement  between  the  two  sets 
of  results.    It  is  likely  that  the  agreement  can  be  further  improved  by  increasing 
the  number  of  panel  from  100,  used  for  these  calculations,  to  200  or  200  more. 
As  was  done  previously  for  an  airfoil  oscillating  in  plunge,  we  will    now 
consider  the  case  of  an  airfoil  oscillating  in  pitch  with  particular  emphasis  on 
drag  or  propulsive  force.    We  will  first  look  at  the  case  without  presence  of 
layers  of  adjacent  wake,  then  look  at  the  wake  effect  including  the  important 
consideration  of  phasing.  We  will  begin  by  considering  the  numerical  results  and 
comparing  these  results  with  Garrick,  Bosch,  for  "no  wake"  case,   and  then  com- 
pare   with  Loewy' s  theory     when     the  wakes  are  included. 
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Considered  is  the  case  of  an  NACA0007  airfoil  set  at  zero  degree  angle 
of  attack.  The  airfoil  is  oscillating  with  an  amplitude    1    degree.  The  reduced 
frequency    is  taken  at  k  =  0.0617  and    the  non-dimensional  wake  spacing  was 
h  =200.  which  is  numerically  equivalent  to  having  the  wake  at  infinity-.  The  time 
history  of  the  lift  variation,  CI,  and  the  time  history  of  drag  variation.  Cd,  reveal 
the  following.    We  observe  that  the  lift  varies  as  the  frequency  of  oscillation, 
with  a  lift  mean  value  of  zero  as  we  expect. 

In  contrast  to  the  lift  time  history,  the  drag  time  history  indicates  that  the 
calculated  drag  vanes  at  twice  the  frequency  of  the  pitch  oscillation.  The  am- 
plitude is  given  by  Cd  =  0.000192.  The  mean  drag  value  is  found  to  be  Cd  = 
0.00009  and  the  mean  lift  is  of  course  0.0  at  zero  angle  of  attack. 

For  an  airfoil  oscillating  in  pitch  with  the  wake  layer  near  the  reference 
airfoil  at  non-dimensional  wake  spacing,  h  =  2.  The  case  of  zero  degree  phase, 
where  the  wake  below  the  airfoil  at  h  =  2  is  in  phase  the  wake  shed  from  the 
reference  airfoil  itself.  Review  of  the  results  indicates  a  reduction  in  the  ampli- 
tude of  lift  variation  for  the  m  =  0.0  case  when  compared  to  the  case  with  no 
wake  at  all  (wake  at  infinity).  We  see  that  the  lift  amplitude  decreases.  We  also 
see  an  increase  in  mean  steady  drag  accompanies  the  reduction  in  lift  amplitude, 
for  the  case  of  zero  wake  spacing  (m  =  0).  Proper  phasing  is  achieved  at  m 
=0.5.  At  this  value,  we  found  that  the  lift  enhancement  in  this  case  of  pure 
pitch,    is    similar  to    what  we  found  in    the    plunge  case. 
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B  CONCLUSIONS 

The  most    important  conclusions  that  can  be  drawn  from  this  research  are; 

1.  For  drag,  in  both  cases  of  pure  plunge  and  pure  pitch,  wake  phasing  can 
increase    or  decrease  the  steady    component  of  drag  acting  on  the  airfoil.     The 
increase  in  the  steady  component  of  propulsive  force  due  to  wake  phasing  in  the 
case  of  pure  plunge  was  sufficient  to  significantly  enhance  the  "  Katzmayr  " 
effect  value  of  propulsive  force. 

2.  For  an  airfoil  oscillating  in  plunge  or  pitch,  for  optimum  reduction  in  drag 
and  enhancement  in  vibratory  lift  the  phase  angle  of  the  wake  vorticity  in  the 
single  wake  case  with  respect  to  motion  of  the  reference  airfoil  should  be  about 
180° 

3.  For  an  airfoil  oscillating  in  plunge  or  pitch,  the  largest  value  of  steady  drag 
and  smallest  value  of  oscillatory  lift  occurs  when  the  phase  angle  of  the  wake 
vorticity    with  respect  to  motion  of  the  reference  airfoil  is  about    0  °. 

4.  For  lift,  in  both  plunge  and   pitch  in  the  single  wake  case,  wake  phasing  can 
increase  or  reduce  the  oscillatory  lift  acting  on  the  airfoil.  Where  in  the  absence 
of  wake  vorticity,  the  effect  of  oscillations  decreases  the  lift  (lift  deficiency  ). 
With    the  wake  present  we  also    observe  lift    enhancement     (  lift  efficiency  ). 

5.  There  exists  a  phase  angle,  75°  for  the  case,  (reduced  k  =  0617,  wake 
spacing,  h  =  2  ),  where  the  effect  of  wake  vorticity  on  the  reference  airfoil  is 
identical  to  the  case  of  the  wake  at  infinity  (  represented  in  the  code  by  wake 
h  =  200  ),  or  no  wake  vorticity  at  all. 
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Applying  these  findings  to  the  OH-6A  case,    we  can  add  the  following 
conclusions  : 

6.  Measured  reductions  in  power  are  feasible  when   the  "Katzmayr  "  effect   and 
thesis  -  reported    wake  enhancements  are  included. 

7.  Measured    higher  harmonic  control  results  show  least  benefit  at  100  knots 
due  to  the  wake  being  transported  greatest  distance  from  rotor  at  this  speed. 

8.   Poor  repeatability  of  measured  open  loop  performance  data  is  to  be  expected 
due  to  shifts  in  wake  position  ! 

9.      For  helicopters,  it  is  standard    rotor    blade  design  to  locate  the  center  of 
gravity,  aerodynamic  center,  and  elastic  axis  at  the  quarter  chord  point  to  avoid 
flutter.  This  was  illustrated  in  the  analysis.    Results  of  this  thesis  showed  that 
flutter  is    impossible  for  this    condition  without  wake  interference.     However, 
Loewy's  results  introduce  the  possibility  of  wake  induced  flutter.  Future  investi- 
gations, therefore,  will  have  to  explore  this  possibility  in  more  detail. 


E.  R.  Wood.    Max  F.  Platzer.    Ahmed  Ahourahma.    Mark    A.  Couch 
On  the  Unsteady  Aerodynamics    of   Higher  Harmonic  Control. 
Paper  No.  C  17.  Nineteenth  European  Rotorcratt  Forum.  Cernobbio 
(Comoi.  Italy.  September  1993. 
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C    RECOMMENDATIONS 

1.  A  new    flight  test  program,  directed    at    studying  helicopter  performance,  is 
needed.  The  NASA-Army    OH-6A  program  was  dedicated  to  vibration  testing. 

2.  Continued  full  scale  HHC  tests  in  the  NASA  Ames  40  '  x  80'  wind  tunnel  is 
encouraged.  In  addition,  it  is  suggested  that  tests  be  initiated  in  small  university 
wind  tunnels  to  verify    the  findings  of  this  thesis. 

3.  Controlled  whirl  tower  testing  is  also  recommended,  since  this  allows  careful 
measurement  and  control  of  performance  parameters. 

4.  A  wind  tunnel  program  is  needed  to  explore  similar  cases  and  to  verify  the 
present  research  findings. 

5.  The  panel  method  is  a  significant  analysis  tool  that  can  be  applied  to  study  and 
solve  similar  problems.  It  provides  an  excellent  method  for  studying  single  or 
two-degree-of-freedom  flutter  problems. 
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